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ABSTRACT 


Three-Dimensional  Response  of  an  Automobile  to 
a  Generalized  Impulse.  (December  1979) 

John  Frederick  Betz,  B.S.,  United  States  Air  Force  Academy 
Chairman  of  Advisory  Cormiittee:  Dr.  Hayes  E.  Ross,  Jr. 

Roadside  design  standards  have  not  kept  pace  with  the  recent  up¬ 
swing  of  compact  and  subcompact  vehicle  traffic  on  the  highways. 
Written  when  the  predominant  vehicle  size  was  larger  than  the  present 
trends,  these  standards  require  examination  in  light  of  the  smaller, 
lighter  vehicles.  By  modifying  the  highway-vehicle-object  simulation 
model  (HVOSM)  ,  a  three-dimensional  computer  simulation  cf  signpost 
impacts  is  available  as  a  tool  to  evaluate  the  standards. 

In  the  course  of  this  thesis,  validation  of  the  modified  HVOSM 
was  completed  to  prove  the  reliability  of  the  program  in  simulating 
real  collisions;  then,  the  parameter  studies  were  carried  out.  Sign¬ 
post  models  meeting  the  upper  limits  of  the  AASHTO  standards  iX&j  were 
prepared.  The  two  models  represented  the  major  signpost  installation 
methods:  a  breakaway  signpost  connected  to  a  base  at  ground  level, 
and  a  base-bending  signpost  placed  directly  into  the  ground.  Then, 
large,  medium,  compact  and  subcompact  vehicle  impacts  with  the  sign¬ 


post  models  were  simulated  using  the  modified  HVOSM 


t; 


In  each  case,  the  subcompact  vehicle  experiences  extreme  accel 


erations  and  total  loss  of  control.  The  breakaway  sign  model  causes 
spinout  of  the  subcompact  car,  while  the  base-bending  model  impact 
results  in  vehicle  rollover.  Each  study  reveals  the  anticipated 


iv 


results:  the  occupants  of  the  two  larger  vehicles  survive  impact  in  a 
stable,  safe  fashion,  while  the  two  smaller  cars,  especially  the  sub¬ 
compact,  experience  hazardous  accelerations  and  vehicle  instability 
after  Impact.  These  results  all  point  to  a  re-evaluation  of  roadside 
design  standards  as  presented  by  AASHTO. 
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CHAPTER  I 

PROBLEM  BACKGROUND 

Introduction 

A  quick  inspection  of  current  attitudes  and  changing  lifestyles  is 
very  revealing  about  the  trend  in  size  today.  The  concept  of  "the  boys 
want  less"  is  even  carried  into  the  automobile  market,  where  the  rising 
prices  and  concurrent  shortages  of  gasoline  have  combined  to  make  the 
large,  old-fashioned  status  symbol  obsolete.  Cars  labeled  as  "gas- 
guzzlers"  are  avoided,  and  the  newer  cars  are  much  smaller  and  lighter 
to  meet  the  increasing  demand  for  higher  range  and  gas  mileage. 

The  national  highway  standards,  however,  do  not  reflect  these 
changes.  Although  any  changes  in  national  standards  should  not  be 
made  to  reflect  mere  trends,  a  quick  examination  of  Europe,  where  gaso¬ 
line  has  been  high-priced  for  some  time,  indicates  that  the  smaller 
cars  will  become  a  larger  proportion  of  traffic  in  the  future. 

The  signposts  and  light  poles  along  the  highways  in  the  U.S.  have 
been  designed  using  existing  AASHTO  specifications  (10),  which  were 
developed  for  compact  sized  (2250  pounds)  automobiles.  The  present 
car  size  trend,  however,  includes  many  subcompact  cars  weighing  1800 
pounds  or  less.  This  is  one  factor  which  indicates  the  need  for  a 
re-examination  of  the  AASHTO  standards. 

Testing  performed  by  the  Texas  Transportation  Institute  (TTI)  (7) 
also  presented  some  interesting  information  to  substantiate  the  need 

The  model  journal  chosen  for  this  thesis  was  the  Journal  of  the 
Engineering  Mechanics  Division,  Proceedings  of  the  American  Society  of 
Civil  Engineers. 
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for  re-evaluation.  These  tests  were  originally  performed  to  test  the 
effects  of  different  signpost  configurations  on  Chevrolet  Vegas,  a 
compact  car,  in  simulated  impacts.  An  interesting  development  in  the 
tests  was  that  the  Vega  came  quite  close  to  rollover  in  several  of  the 
impacts.  The  vehicle  actually  did  experience  rollover  in  one  test 
(test  #11),  but  the  sign  configuration  did  not  meet  existing  AASHTO 
standards.  However,  the  important  fact  was  that  the  rollover  was 
nearly  the  result  in  several  collisions  with  signs  that  met  the  stan¬ 
dards.  Smaller  cars  are  more  likely  to  experience  rollover  than  larger 
cars,  so  it  is  quite  possible  that  subcompact  cars  would  have  actually 
rolled  over  in  the  same  situations.  This  second  factor,  the  near  roll¬ 
over  of  compact  cars  in  impacts  with  signposts  meeting  AASHTO  stan¬ 
dards,  combined  with  the  car  size  trend  toward  smaller  and  lighter 
automobiles,  formed  the  basis  for  this  thesis. 

The  concern  of  this  thesis  is  then  the  examination  of  signpost  im¬ 
pacts  using  computer  simulation.  The  behavior  of  subcompact,  compact, 
medium,  and  large  cars  in  collisions  with  a  signpost  or  light  pole 
meeting  AASHTO  design  criteria  will  be  studied  to  reveal  whether  the 
existing  standards  allow  rollover  to  occur  in  the  smaller  automobiles. 

Program  Selection 

Two  existing  computer  simulations  of  automobiles  were  examined  to 
determine  which  was  better  suited  to  study  signpost  impacts.  These 
programs  were  the  Highway-Vehicle-Object  Simulation  Model  (HVOSM)  (8) 
and  the  Simulation  Model  of  Automobile  Collisions  (SMAC)  (4).  Each  one 
had  negative  and  postive  aspects  which  were  examined  in  light  of  the 
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requirements  of  the  planned  simulations.  The  final  choice  turned  out 
to  be  much  easier  than  anticipated  due  to  the  strong  negative  aspects 
of  one  of  the  programs. 

The  HVOSM  was  developed  to  simulate  the  three-dimensional  motion 
of  an  automobile  subjected  to  roadside  guardrail  collisions  or  a  run¬ 
off-road  accident.  It  was  also  developed  to  study  the  ride  and 
handling  behavior  of  an  automobile  traversing  a  given  roadway  configu¬ 
ration.  The  program  has  been  effectively  used  in  guardrail  evalua¬ 
tions  and  roadside  design  for  several  years.  The  shortcoming  of  the 
current  program  is  in  the  lack  of  provisions  to  simulate  impacts  of 
the  model  vehicle  with  roadway  signs  or  lampposts.  Thus,  if  this 
program  were  to  be  used  in  this  study,  extensive  modifications  would 
be  required  to  add  impulse  evaluation  to  the  possible  program  func¬ 
tions.  However,  the  program  has  the  advantages  of  three-dimensional 
motion  and  proven  performance  in  the  design  of  highways  and  roadside 
structures. 

SMAC  is  another  computer  program,  and  it  is  used  to  evaluate  im¬ 
pacts  between  automobiles  and  signposts.  SMAC  was  originally  developed 
by  Calspan  and  further  modified  by  TTI.  Despite  the  fact  that  the 
program  was  developed  to  study  signpost  impacts,  there  are  two  inherent 
weaknesses  in  the  program  which  invalidate  its  use  in  this  study.  First 
of  all,  the  program  only  considers  two-dimensional  (planar)  motion. 

Since  this  study  is  concerned  with  possible  vehicle  rollover,  a  planar 
model  of  collision  is  unsuitable.  Finally,  the  program  assumes  that 
the  signposts  are  perfectly  rigid.  Breakaway  posts  are  designed  to 
prevent  the  negative  aspects  of  rigidity,  and  nonbreakaway  posts  create 
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large  downward  forces  on  automobiles  when  the  signposts  plastically 
deform  and  wrap  around  the  car  frame  as  the  vehicle  rides  them  down. 
Since  the  SMAC  model  does  not  consider  the  three-dimensional  motion 
required  for  rollover  and  it  evaluates  only  rigid  posts  in  collisions, 
the  SMAC  model  is  unsuited  for  this  study.  The  modification  required 
for  the  SMAC  model  would  amount  to  changing  the  entire  program,  which 
would  be  a  greater  task  than  modifying  the  HVOSM. 

Decision 

An  analysis  of  the  features  of  the  HVOSM  and  the  SMAC  programs 
resulted  in  the  choice  of  the  HVOSM  as  the  program  best  suited  for 
this  study.  Modifications  had  to  be  made  to  the  HVOSM  to  allow  the 
option  of  accepting  generalized  impulse  data  to  simulate  signpost  im¬ 
pacts.  The  program  was  to  simulate  the  effects  of  any  general  force 
acting  on  any  location  of  the  vehicle.  Since  any  general  force  can 
be  broken  down  into  its  components  aligned  with  the  vehicle  axes, 
the  program  would  use  the  time-dependent  force  and  location  data  as 
referenced  to  and  oriented  with  the  principle  vehicle  axes.  The 
location  data  would  locate  the  point  of  the  force  application  from 
the  sprung  mass  center  of  gravity  of  the  vehicle  measured  along  those 
axes.  The  implementation  of  these  modifications  to  the  HVOSM  is  the 
subject  of  the  next  chapter.  Further  references  and  a  figure  relating 
to  the  HVOSM  can  be  found  in  Appendix  II. 
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CHAPTER  II 

PROGRAM  MODIFICATION 
Background 

As  previously  mentioned,  the  HVOSM  is  a  proven  program  of  numerous 
capabilities.  Mainly  used  for  roadside  barrier  design  and  studies  of 
automobile  behavior  on  varying  surfaces,  the  program  had  to  be  modified 
to  accept  new  options  without  interfering  with  its  existing  abilities. 
After  correcting  the  errors  due  to  differences  in  the  compilers  (Appen¬ 
dix  III),  the  program  was  studied  at  length  to  insure  that  the  changes 
were  made  to  affect  the  fewest  subroutines  and  therefore  minimize  any 
possible  errors.  The  changes  themselves  were  quite  simple,  and  the 
whole  series  of  modifications  was  simply  a  matter  of  finding  where  the 
critical  operations  took  place. 

The  changes,  and  thus  the  divisions  of  this  chapter,  can  be 
classed  into  three  areas:  the  input  phase,  the  mathematical  opera¬ 
tions,  and  the  output  phase.  The  actual  computer  listings  of  the 
changes  as  they  appear  in  the  program  can  be  found  in  Appendix  IV.  The 
theory  and  implementation  of  the  changes  are  the  subject  of  the  divi¬ 
sions  in  this  chapter. 

Input  Phase 

The  HVOSM  had  to  be  modified  to  accept  a  new  matrix  of  time- 
dependent  force  and  location  data.  Fifty  time  points  were  considered 
to  be  adequate  to  describe  a  sign  impact  with  the  use  of  linear  inter¬ 
polation  between  time  points.  Since  the  matrix  had  to  define  the 
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time,  the  forces  in  the  three  principal  directions,  and  the  three- 
dimensional  location  of  the  point  of  application,  the  matrix  was  sized 
as  7  by  50.  Another  necessary  pair  of  additions  were  a  variable  that 
could  be  used  to  trigger  the  sign  impact  program  option  into  opera¬ 
tion,  and  a  variable  to  define  the  total  number  of  time  points  used  for 
each  problem  (since  fifty  time  points  would  not  always  be  necessary). 

As  each  of  these  new  variables  were  to  be  used  in  a  number  of  sub¬ 
routines,  they  were  added  to  existing  common  blocks  in  the  program 
rather  than  adding  to  the  argument  of  the  individual  subroutine.  The 
variables  were  named  as  follows:  The  time,  force,  and  location  matrix 
was  named  SUBB(7,50);  the  sign  impact  option  trigger  was  named  COLL; 
and  the  variable  for  the  number  of  time  points  was  named  NUMT  (an 
integer- valued  variable).  The  SUBB  matrix  and  the  COLL  variable  were 
then  added  to  the  end  of  the  BARIER  common  block  throughout  the  pro¬ 
gram,  and  the  NUMT  variable  was  added  to  the  INPT1  common  block 
throughout  the  program.  The  initial  values  of  these  values  were  all 
cleared  to  zero  in  the  initial  part  of  the  main  program.  COLL  and  the 
SUBB  matrix  were  cleared  along  with  the  rest  of  the  BARIER  common 
block  using  the  ERASE  subroutine,  while  NUMT  was  cleared  to  zero  by 
adding  a  simple  program  statement. 

After  the  variables  were  set  up  and  ready  to  accept  the  proper 
numbers,  the  actual  input  operation  was  next  completed.  This  involved 
modifying  the  INPUT  subroutine  to  accept  the  additional  values.  In 
the  HVOSM,  the  input  is  organized  according  to  ICARD  values  which  are 
specified  by  a  number  in  columns  79  and  80  on  the  input  data  cards. 

The  ICARD  values  are  used  to  specify  a  GO  TO  statement  location  in  the 
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INPUT  subroutine,  where  the  values  from  the  cards  are  affixed  to  the 
proper  variables.  The  cards  are  divided  into  fields  of  eight  charac¬ 
ters  each,  with  nine  fields  per  card.  Thus,  some  ICARD  values  did  not 
use  the  full  allotment  of  nine  fields,  and  these  unused  fields  could 
be  used  to  define  the  COLL  and  NUMT  variable  values.  The  great  size  of 
SUBB  matrix  warranted  creation  of  a  new  ICARD  series  to  transfer  the 
matrix  values. 

A  further  explanation  of  the  workings  of  the  INPUT  subroutine  will 
help  to  ease  the  explanation  of  the  changes  performed  on  the  subroutine 
to  load  these  variables.  The  subroutine  reads  the  data  cards  and  re¬ 
cords  them  in  card  images  which  are  then  printed  in  the  output  phase 
so  that  the  values  as  the  computer  reads  them  can  be  checked  for  error. 
However,  those  same  values,  which  are  contained  in  nine  fields  of  eight 
characters  each,  are  loaded  into  nine  locations  in  a  temporary  one¬ 
dimensional  matrix  named  DUM(9).  The  values  from  each  of  the  nine 
fields  are  loaded  into  DUM(l)  through  DUM(9),  respectively.  Then,  the 
GO  TO  statement,  keyed  by  the  ICARD  value,  leads  the  program  to  the 
correct  section,  where  those  values  are  loaded  from  the  DUM  matrix 
into  the  corresponding  program  variables.  If  an  ICARD  series  does  not 
use  the  full  nine  fields,  then  some  of  the  DUM  matrix  locations  are 
not  utilized  and  these  openings  are  available  for  variable  additions 
without  changing  the  program  (other  than  the  required  variable  loading 
statements  with  the  unused  matrix  locations). 

The  ICARD=2  data  set  only  used  seven  of  the  nine  fields  on  the 
data  card,  so  the  eighth  and  ninth  fields  on  that  card  were  used  to 
define  COLL  and  NUMT,  respectively.  Then,  since  the  ICARD=2  value 


sends  control  to  statement  number  4  in  the  INPUT  subroutine,  the  state¬ 
ments 

COLL  =  DUM(8) 

NUOT  =  DUM(9) 

were  added  after  the  statement  assigning  the  DUM(7)  value  to  a  variable 
name.  Thus,  simply  by  using  the  full  capabilities  of  the  ICARD=2  data 
set,  COLL  and  NUMT  input  values  were  added  to  the  program  with  only  two 
program  statements. 

The  SUBB  matrix  presented  a  greater  challenge.  First  of  all,  each 
of  the  matrix  locations  had  to  be  defined  for  the  proper  input  align¬ 
ment.  The  locations  of  the  SUBB(7,50)  matrix  were  defined  as  follows, 
using  the  general  location  SUBB(A,B): 

1.  All  A=1  locations  are  discrete  time  values,  in  seconds. 

2.  All  A=2  locations  are  instantaneous  force  magnitudes  aligned 
with  the  local  vehicle  x-axis,  in  pounds  (+  or  -). 

3.  All  A-3  locations  are  instantaneous  force  magnitudes  aligned 
with  the  local  vehicle  y-axis,  in  pounds  (+  or  -). 

4.  All  A=4  locations  are  instantaneous  force  magnitudes  aligned 
with  the  local  vehicle  z-axis,  in  pounds  (+  or  -). 

5.  All  A=5  locations  are  the  distance  from  the  sprung  mass  center 
of  gravity  to  the  point  of  force  application  measured  along 
the  vehicle  x-axis,  in  inches  (+  or  -). 

6.  All  A=6  locations  are  the  distance  from  the  sprung  mass  center 
of  gravity  to  the  point  of  force  application  measured  along 
the  vehicle  y-axis,  in  inches  (+  or  -). 

7.  All  A=7  locations  are  the  distance  from  the  sprung  mass  center 
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of  gravity  to  the  point  of  force  application  measured  along 
the  vehicle  z-axis,  in  inches  (+  or  -). 

8.  All  B  values,  from  a  minimum  of  2  up  to  a  possible  50,  are  the 
values  of  each  set  of  forces  and  locations  corresponding  to 
each  discrete  time  value. 

As  can  be  seen  from  the  definitions  listed  above,  the  forces  all 
act  at  a  point  for  each  time  value. 

Next,  the  ICARD=29  provision  (the  last  program  data  set  ICARO 
value  had  been  28)  had  to  be  added  to  the  GO  TO  statement,  using  the 
statement  number  chosen  for  the  additional  data  set.  Statement  number 
53  was  not  in  use  elsewhere  in  the  INPUT  subroutine,  so  it  was  used  to 
locate  the  ICARD=29  variable  loading  statements.  After  adding  location 
53  to  the  GO  TO  argument,  statement  number  53  and  the  rest  of  the  ICARD 
=29  variable  loading  statements  were  inserted  directly  between  the 
next-to-last  statement  and  the  END  statement  in  the  INPUT  subroutine. 

All  that  remained  was  to  add  the  statements  to  transfer  the  values 
from  the  DUM  matrix  positions  to  the  proper  locations  in  the  SUBB  ma¬ 
trix.  An  additional  counter  variable  was  added  for  this  purpose.  The 
variable  JJ  was  set  equal  to  zero  at  the  start  of  the  program  with  an 
additional  statement.  JJ  defines  the  location  of  each  set  of  data  for 
each  time  set,  thereby  sequentially  loading  the  seven  variables  into 
the  SUBB  matrix.  This  requires  that  the  ICARD=29  data  deck  be  orga¬ 
nized  in  order  of  increasing  time  values,  both  for  this  sequential 
loading  and  for  use  in  linear  interpolation  later  in  the  program.  It 
is  vitally  important  to  the  operation  of  the  program  that  the  data 
cards  be  in  order  of  sequentially  increasing  time  values. 
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After  defining  the  SUBB  matrix  and  inserting  JJ  as  a  counter 
variable,  the  following  statements  were  added  to  complete  the  input  of 
the  SUBB  matrix: 

53  JJ  =  JJ+1 

SUBB(l.JJ)  =  DUM(l) 

SUBB(2,JJ)  =  DUM(2) 

SUBB(3,JJ)  =  DUM(3) 

SUBB(4,JJ)  =  DUM(4) 

SUBB(5,JJ)  =  DUM ( 5 ) 

SUBB(6,JJ)  =  DUM(6) 

SUBB(7  ,JJ)  =  DUM(7) 

GO  TO  2 

Therefore,  the  data  card  only  uses  seven  of  the  nine  fields  provided, 
and  the  variables  are  loaded  in  the  order  of  time,  x,  y,  and  z  forces, 
and  the  x,  y,  and  z  locations. 

That  was  the  extent  of  the  necessary  modifications  to  input  the 
values  for  the  sign  impact  option.  The  important  changes  to  remember 
to  operate  the  sign  impact  option  are  that  COLL  and  NUMT  are  the  eighth 
and  ninth  variables  on  the  ICARD=2  data  set,  where  COLL  is  input  as  1.0 
to  trigger  the  sign  impact  option  and  NUMT  is  set  equal  to  the  number 
of  discrete  time  points  used  in  the  SUBB  matrix  to  define  the  impulse 
input.  Also,  the  ICARD=29  data  set,  which  must  be  entered  in  the 
order  of  sequentially  increasing  time  steps,  is  used  to  load  the  SUBB 
matrix.  The  first  seven  fields  are  used,  defined  as  the  time  (in 
seconds),  the  x,  y,  and  z  forces  (in  pounds),  and  the  x,  y,  and  z  lo¬ 
cations  (in  inches),  respectively,  in  the  first  through  seventh  fields. 
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Important  requirements  are  that  the  first  time  step  should  always  be 
zero,  and  the  last  time  step  should  be  set  after  the  program  is  antici¬ 
pated  to  stop  running.  In  other  words,  if  the  program  is  set  to  inves¬ 
tigate  only  the  first  3  seconds  of  car  motion,  then  the  last  time  step 
should  be  set  to  greater  than  3  seconds  for  proper  program  operation. 

Mathematical  Operations 

After  the  time,  force,  and  location  data  and  the  trigger  variable 
were  input  to  the  program,  the  next  obstacle  was  to  find  exactly  where 
in  the  program  to  insert  the  equations  used  to  compute  the  forces  and 
moments  acting  on  the  automobile  during  sign  impact.  A  thorough  analy¬ 
sis  of  the  program  provided  the  location  which  required  relatively  few 
additions  to  the  program  to  add  the  necessary  equations. 

The  SFORCE  subroutine,  with  the  help  of  several  subordinate  sub¬ 
routines,  calculates  the  forces  and  moments  acting  on  the  vehicle  due 
to  roadside  guardrail  collisions.  SFORCE  calls  several  subordinate 
subroutines  in  the  course  of  its  operation,  most  notably  RESFRC,  which 
uses  the  information  developed  in  the  SFORCE  subroutine  to  compute  the 
forces  and  moments  acting  on  the  automobile. 

Since  all  these  operations  could  be  narrowed  down  to  the  SFORCE 
subroutine,  the  sign  impact  modifications  were  placed  at  the  tail  end 
of  the  subroutine,  and  the  COLL  trigger  was  placed  to  circumvent  the 
normal  operation  of  SFORCE  and  the  subordinate  subroutine.  As  all 
the  collision  forces  and  moments  used  by  the  HVOSM  originated  in  the 
SFORCE  subroutine  and  its  subordinate  subroutines,  the  addition  of  the 
signpost  impact  option  trigger  and  the  correct  equations,  using  the 
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same  variable  names  used  by  the  RESFRC  subroutine,  to  the  SFORCE  sub¬ 
routine  were  the  only  necessary  program  modifications  for  the  mathema¬ 
tical  operations. 

First  of  all,  the  SFORCE  subroutine  sign  impact  option  trigger 
statement  had  to  be  set  up  and  placed  before  the  first  operational 
subroutine  statements.  The  COLL  variable  was  used  in  its  trigger  func¬ 
tion  at  the  head  of  the  program,  directly  following  the  clearing  of  the 
variables  (when  they  are  set  equal  to  zero).  The  statement  added  was 

IF(C0LL.EQ. 1. )  GO  TO  1002 

where  1002  is  the  statement  number  of  the  first  of  the  sign  impact 
modification  option  statements.  As  defined  in  the  input  phase,  COLL  is 
set  equal  to  1.0  to  trigger  the  sign  impact  option  in  the  ICARD=2  data 
set. 

Finally,  the  statements  for  the  mathematical  operations  had  to  be 
organized  and  placed  starting  at  statement  number  1002.  Also,  the  pro¬ 
gram  variables  for  the  forces  and  moments,  oriented  with  the  local 
vehicle  axes  and  in  the  proper  program  units,  had  to  be  used  to  insure 
the  rest  of  the  subroutines  would  smoothly  continue  their  regular  func¬ 
tions.  The  program  units  for  the  forces  were  found  to  be  pounds,  and 
the  moments  were  found  to  be  in  inch-pounds.  By  specifying  the  input 
values  for  the  forces  and  locations  in  the  units  of  pounds  and  inches, 
respectively,  this  requirement  had  already  been  met.  Next,  the  pro¬ 
gram  variable  names  for  the  x,  y,  and  z  forces  and  the  x,  y,  and  z 
moments  in  RESFRC  were  found  to  be  SFXS,  SFYS,  SFZS,  and  SNPS,  SNTS, 
and  SNPSS,  respectively. 

All  that  remained  was  to  set  up  the  linear  interpolation  for  the 
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forces  and  locations  at  time  values  between  the  input  time  values.  The 
program  variable  for  time  was  found  to  be  T.  Using  the  difference  be¬ 
tween  T  and  the  next  smaller  time  step,  the  variable  VART  was  created 
as  a  general  percentage  value  from  the  linear  interpolation  method. 

The  mathematical  operation  proceeded  as  follows:  First  of  all, 
the  T  which  the  program  had  progressed  to  had  to  be  situated  between 
the  time  step  greater  than  T  and  the  next  smallest  time  step.  In  other 
words,  the  location  of  T  in  reference  to  the  input  time  values  in  SUBB 
had  to  be  found.  That  is  why  the  sequential  loadings  of  the  time 
values  had  to  be  performed,  because  the  proper  orientation  of  the  T 
variable  used  by  the  program  could  not  be  performed  unless  the  time 
values  were  in  sequentially  increasing  order.  After  the  location  of  T 
was  found  using  an  IF  statement  loop,  the  next  operation  was  to  compute 
VART  as  a  percentage  value  using  T.  This  was  done  by  dividing  the  dif¬ 
ference  between  T  and  the  next  smallest  time  step  by  the  difference  be¬ 
tween  the  next  largest  time  step  and  the  next  smallest  time  step  (both 
relative  to  T).  Then,  this  percentage  was  used  to  compute  the  forces 
in  the  x,  y,  and  z  directions  at  the  time  T.  Next,  the  x,  y,  and  z 
locations  at  time  T  were  found  using  the  same  linear  interpolation 
method  (and  called  ECKS,  WHY,  and  ZEE,  respectively). 

After  computing  the  values  for  SFXS,  SFYS,  and  SFZS  and  the  loca¬ 
tion  values  ECKS,  WHY,  and  ZEE  at  time  T,  the  operations  were  con¬ 
cluded  by  using  simple  equations  to  compute  the  moments  SNPS,  SNTS, 
and  SNPSS.  The  whole  operation,  in  computer  language,  was  added  as 
f ol 1 ows : 
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1002  I  =  2 

1003  IF(T.LE.SUBB(1,I) )  GO  TO  1004 
I  =  1+1 

GO  TO  1003 

1004  IM1  =  1-1 

VART  =  (T-SUBB(i, IM1 ) ) / ( SUBB (1,1) -SUBB ( 1 , IM1) ) 

SFXS  =  SUBB(2, IM1)+VART*(SUBB(2, I )-SUBB(2 , IM1) ) 

SFYS  =  SUBB(3, IM1)+VART*(SUBB(3, I )-SUBB(3, IM1) ) 

SFZS  =  SUBB (4 , IM1 )+VART*( SUBB (4,1 ) -SUBB (4 , IM1 ) ) 

ECKS  =  SUBB ( 5 , IM1 ) +VART*( SUBB (5,1) -SUBB ( 5 , IM1 ) ) 

WHY  =  SUBB (6 , IM1)+VART*(SUBB(6 , I ) -SUBB(6 , I Ml ) ) 

ZEE  =  SUBB(7, IM1)+VART*(SUBB(7 ,1 )-SUBB(7 , IM1) ) 

SNPS  =  SFZS*WHY-SFYS*ZEE 
SNTS  =  SFXS*ZEE-SFZS*ECKS 
SNPSS  =  SFYS*ECKS-SFXS*WHY 

So,  the  mathematical  operations  were  confined  to  one  subroutine, 
SFORCE,  and  were  added  without  reducing  any  of  the  existing  HVOSM 
capabilities.  The  only  limitation  created  by  this  method  is  that  the 
HVOSM  program  cannot  simulate  a  guardrail  collision  and  a  signpost  im¬ 
pact  in  the  course  of  the  same  computer  run.  That  is  not  considered  a 
limitation,  however,  due  to  the  rarity  cf  the  event  in  the  normal  pro¬ 
gram  use. 


Output  Phase 

By  placing  the  mathematical  operations  in  the  SFORCE  subroutine 


and  using  program  variables  for  forces  and  moments,  no  additional 


output  statements  were  required  for  the  program  since  the  regular  PRINT 
and  OUTPUT  subroutines  list  all  the  necessary  output  for  the  program 
results.  The  only  additional  output  required  by  the  signpost  modifica¬ 
tion  was  the  listing  of  the  values  input  in  the  SUBB  matrix. 

The  input  variables  are  printed  out  under  statements  in  the  IDOUT 
subroutine.  Therefore,  the  output  modifications  had  to  be  performed  on 
the  IDOUT  subroutine,  once  again  using  the  COLL  trigger  associated  with 
the  sign  impact  option.  The  BARIER  common  block  had  to  be  added  to  the 
IDOUT  subroutine  so  that  the  variables  would  be  present  for  output, 
which  had  been  unnecessary  previously  (since  none  of  the  other  BARIER 
variables  were  derived  directly  in  the  input  operation). 

Proceeding  to  the  end  of  the  IDOUT  subroutine,  the  changes  were 
inserted  immediately  before  the  "CALL  DRIVID"  statement,  first,  the 
COLL  trigger  was  added,  which  transferred  program  flow  to  the  "CALL 
DRIVID"  statement  if  COLL  was  not  equal  to  1.0.  Remember,  COLL  was 
set  equal  to  1.0  for  sign  impact  option  operation  of  the  HVOSM.  Then, 
a  title  line  which  delineated  the  information  set  forth  below  it  as 
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"Sign  Impact  Force  and  Location  Data",  with  a  series  of  astericks,  was 
set  up  for  print.  This  was  followed  by  headings  for  each  column  of 
matrix  information:  time,  forces,  and  locations,  all  with  the  attached 
units.  Rather  than  list  these  statements,  which  are  merely  a  series 
of  formats  and  variable  listings,  reference  is  made  to  Appendix  IV, 
which  contains  the  listings  of  all  the  changes  mentioned  in  this 
chapter  exactly  as  they  appear  in  the  program. 
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Conclusion 

While  at  first  seeming  a  complicated  task,  the  actual  modification 
of  the  HVOSM  was  performed  quite  easily  following  proper  program  analy¬ 
sis.  All  of  the  changes  were  made  without  reducing  the  HVOSM  capabili¬ 
ties  existing  before  the  modifications.  All  of  the  changes  were  made 
in  the  simplest  manner  possible,  and  each  of  the  changes  was  accompanied 
by  comnentary  blocks  to  identify  the  changes  and  their  purpose.  Appen¬ 
dix  IV  contains  the  actual  program  listings  of  the  changes  and  the 
comments  which  precede  them. 
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CHAPTER  III 
PROGRAM  VALIDATION 

General  Approach 

With  the  completion  of  the  HVOSM  program  modification,  the  only 
task  required  before  the  actual  sign  impact  parameter  study  was  the 
validation  of  the  modified  program.  Even  though  the  HVOSM  was  a  proven 
program  before  the  changes  were  made,  the  modified  HVOSM  had  to  demon¬ 
strate  operations  which  could  duplicate  real  conditions  as  closely  as 
possible.  To  that  end,  a  series  of  tests  were  arranged  to  prove  that  r 

the  program  could  match  theoretical  results,  and  also  that  it  could 

approximate  the  results  of  actual  crash  tests  using  limited  crash  test  ; 

data. 

The  series  of  four  validation  exercises  were  arranged  in  order  of 
increasing  complexity.  First  of  all,  a  simple  linear  deceleration 
example  was  performed  to  compare  theoretical  results  with  program  re¬ 
sults.  Secondly,  a  LaGrangian  analysis  of  a  sample  vehicle  with  a 
constant  force  aligned  with  the  vehicle  axis  was  used  in  a  comparison 
between  theoretical  and  program  simulated  planar  motion.  After  the 
successful  completion  of  these  two  tests,  the  next  two  exercises  were 
trials  at  duplicating  actual  crash  test  results.  The  crash  tests  in 
question  were  crash  tests  #11  and  #9,  performed  as  described  in 
reference  7.  While  the  only  data  used  in  duplication  of  those  tests 
were  the  accelerometer  readouts  (along  with  all  the  vehicle  parameters), 
the  program  still  closely  approximated  the  actual  results  of  the  tests. 

A  closer  examination  of  each  of  the  exercises  will  reveal  the  function 
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of  each  in  confirming  the  validity  of  the  modified  HVOSM. 

The  test  vehicle  parameters  for  the  theoretical  tests  were  taken 
from  the  HVOSM  User's  Manual  (3).  These  parameters  were  taken  for  a 
standard  sized  car,  and  the  tests  were  run  using  a  frictionless  road 
both  for  the  theoretical  analysis  and  the  program  simulation.  The 
crash  tests  provided  particular  vehicle  parameters  for  the  automobile 
used  in  the  crash  tests,  so  the  input  data  reflected  a  particular  auto¬ 
mobile  make  and  model  (1971  Vega)  and  a  proper  surface  with  a  friction 
factor  of  .5. 


Linear  Deceleration 

The  linear  deceleration  test  was  a  basic  test  of  the  proper  func¬ 
tioning  of  the  modified  HVOSM.  The  program  was  run  after  theoretical 
analysis  showed  a  required  force  of  6829.02  pounds  applied  for  1.5 
seconds  to  the  center  of  gravity  to  stop  the  sample  vehicle  with  a  for¬ 
ward  speed  of  60  miles  per  hour.  The  theoretical  equation  used  for 
this  simple  test  was: 

MaV  =  FAt 

Concentrating  the  impulse  force  at  the  sprung  mass  center  of  gravity  to 
duplicate  the  lumped  mass  assumption  of  the  theoretical  analysis,  the 
HVOSM  results  were  hopefully  to  agree  very  closely  with  the  theoretical 
results.  The  actual  comparison  of  the  results,  shown  as  a  graph  of  the 
automobile  velocity  versus  time  for  both  the  theoretical  results  and 
the  HVOSM  program  output,  appear  in  Figure  1  as  a  simple  line. 

Figure  1  confirms  the  basic  soundness  of  the  modified  HVOSM.  The 


results  agree  almost  exactly,  with  the  greatest  variation  between 
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program  results  and  the  theoretical  results  only  a  negligible  .02  feet 
per  second.  This  is  why  the  results  of  both  are  represented  by  a 
single  line  —  the  differences  were  too  small  to  appear  on  the  graph. 
The  very  fact  that  the  theoretical  lumped  mass  comparison  agreed  so 
closely  with  the  HVOSM's  complex  automobile  model  is  an  affirmation 
that  the  program  is  basically  sound.  The  subsequent  tests  will  deter¬ 
mine  if  the  more  complex  problems  can  be  handled  by  the  modified  HVOSM 
in  like  manner. 


LaGrangian  Analysis 

With  the  basic  program  functioning  assured,  the  next  test  was  to 
prove  the  program  could  operate  in  two  dimensions  and  approximate  theo¬ 
retical  results.  The  theoretical  problem  for  this  test  was  chosen  as 
a  simple  application  of  a  constant  force  to  a  moving  automobile.  The 
constant  200  pound  force  was  to  be  applied  to  the  front  of  the  automo¬ 
bile  in  the  negative  vehicle  x-axis  direction,  that  is,  aligned  con¬ 
stantly  towards  the  rear  of  the  vehicle.  A  LaGrangian  analysis  of  the 
system  was  then  performed  using  the  test  problem.  The  global  and  local 
coordinate  systems  used  in  the  analysis  are  shown  in  Figure  2. 

LaGrange's  equation,  simply  stated,  is 

A.  (JA\  _9U_  +  _£V_  iSl  n  \ 

3qj  ”  3c!i  3qi  " 

where 

LI  =  kinetic  energy 
V  =  potential  energy 

Qe  =  work  done  by  applied  external  forces 


22 


qj  =  a  generalized  (global)  coordinate 

q-j  =  the  first  time  derivative  of  the  generalized  coordinate 
Expressing  the  equation  variables  in  terms  of  global  coordinates. 


U  =  mv2  +  y  I  e2 


(2) 


v2  =  x'2  +  y'2  (3) 

There  are  no  potential  energy  terms  for  this  problem,  so  that  term 
can  be  excluded  from  the  equation.  That  leaves  the  external  work  term 
to  be  defined  in  terms  of  the  generalized  coordinates: 


VfxxVf/f 

(4) 

r' 

where 

* 

x'p  =  x'  +  Xp  cose  +  yp  sine 

(5) 

! ; 

y'p  =  y'  +  Xp  sine  -  yF  cose 

(6) 

and 

Fx  =  F(t)  cose 

(7) 

Fy  =  F(t)  sine 

(8) 

Now  that  all  the  terms  are  defined,  each  of  the 

three  equations 

can  be  extracted  by  allowing  q.  to  be  equal  to  x',  y' 

LaGrange's  equation  —  that  is,  equation  1. 

,  and  e  in 

U  =  1  m(x'2  +  y'2)  +  |'le2 

(9) 
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(10) 
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d  ,  all 
^9x' 


■)  =  ^{mx1 )  =  mx' 


(13) 


’  &<*'  >  =  04) 

JE<f  ■  ^  ■  «  <’5) 


Combining  equations  4,  5,  and  6  for  the  expression  for  Qg, 

3Qe 

g3jr  =  Fx  =  F(t)  cose  (16) 

3Qp 

=  Fy  =  F(t)  sine  (17) 

3Qe 

—  =  fx(-xF  sine  +  yF  cose)  +  Fy(xp  cose  +  yf  sine)  (18) 

Therefore,  the  three  equations  for  theoretical  analysis  end  up  as 

mx'  =  F(t)  cose  (19) 

my'  »  F(t)  sine  (20) 

10  -  2F(t)yF  (21) 


Since  F(t)  is  constant,  the  e  acceleration  is  constant  and  the  x' 
and  y'  accelerations  are  only  dependent  on  e.  The  rest  of  the  theore¬ 
tical  analysis,  which  consisted  of  vehicle  location,  velocity,  and  yaw 
speed  calculations  at  .02  second  time  intervals,  was  performed  using 
simple  equations  of  motion: 


OLD  +  xAt 

(22) 

OLD  +  x0lDAt  +  1  XAt2 

(23) 

These  equations  were  applied  to  each  of  the  three  variables  in  the  cor¬ 
rect  order,  since  the  accelerations  in  the  x’  and  y'  directions  were 
dependent  on  the  value  of  9  from  the  previous  time  step. 


\ 
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After  the  LaGrangian  analysis  for  the  constant  200  pound  force  was 
completed  using  the  same  vehicle  parameters  that  were  used  for  the 
HVOSM  program  input,  the  results  of  each  form  of  analysis  were  graphi¬ 
cally  compared  in  Figures  3  through  6.  Each  of  the  graphs  differ  only 
after  a  significant  period  of  time,  and  even  then  the  differences  are 
very  small.  Figure  3,  the  comparison  of  velocities  in  the  local  vehi¬ 
cle  x-direction,  shows  the  results  differing  by  1.57  feet  per  second 
at  the  end  of  1.5  seconds,  an  error  of  only  3.8%.  The  graph  of  the 
velocities  in  the  vehicle  y-direction,  shown  in  Figure  4,  reveals  a 
maximum  difference  of  1.31  feet  per  second  between  the  two  methods  of 
analysis,  an  error  of  2%.  In  Figure  5,  which  is  the  graph  of  the 
vehicle  location  along  the  global  x-axis,  the  results  agreed  so  closely 
that  only  one  line  could  be  presented.  The  maximum  difference  was  only 
1.33  inches  in  214.8  feet,  representing  virtually  exact  agreement. 
Finally,  Figure  6,  the  vehicle  yaw  speed,  shows  a  maximum  difference  of 
1.18  degrees  per  second,  an  error  of  only  2.6%. 

These  results  served  to  complete  the  validation  of  the  modified 
HVOSM  in  reference  to  theoretical  analysis.  Using  linear  deceleration 
and  LaGrangian  analysis,  the  program  was  examined  for  soundness  in  both 
linear  and  planar  motion.  The  LaGrangian  analysis,  involving  two 
dimensions,  provided  a  number  of  variables  which  were  used  to  evaluate 
the  HVOSM  program  results.  Each  comparison  served  to  further  support 
the  validity  of  the  modified  HVOSM,  revealing  only  negligible  error 
attributed  mainly  to  the  lumped  mass  approximations  used  for  the  iner¬ 
tial  properties  in  the  theoretical  anlaysis. 

Since  the  HVOSM  was  developed  from  theory,  reliable  duplication  of 
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theoretical  results  was  expected.  The  sprung  mass,  with  its  six  degrees 
of  freedom,  and  the  vertical  deflections  of  the  four  unsprung  mass 
(wheels)  provide  only  ten  degrees  of  freedom,  with  a  possible  eleventh 
degree  of  freedom  using  the  free  steering  mode  of  program  operation  (as 
all  these  tests  did).  The  next,  more  difficult  hurdle  was  the  duplica¬ 
tion  of  results  of  actual  crash  tests,  with  their  three-dimensional 
motion  and  infinite  degrees  of  freedom. 

Crash  Test  Number  11 

Duplication  of  real-world  signpost  collisions  is  the  central  func¬ 
tion  of  the  modified  HVOSM.  Necessarily,  the  program  had  to  prove  it¬ 
self  by  duplicating  the  results  of  previous  crash  tests  before  it  could 
be  used  with  confidence.  The  first  crash  test  chosen  for  this  final 
phase  of  validation  was  crash  test  #11  (7)  as  performed  by  TTI .  Accel¬ 
erometers  had  been  placed  on  the  vehicle  and  those  readouts  were  avail¬ 
able,  along  with  the  change  in  momentum  at  .183  seconds  after  the 
collision  and  the  approximate  time  that  vehicle  rollover  started  (this 
vehicle  rolled  over  in  the  test).  The  fact  that  the  vehicle  actually 
rolled  over  was  an  important  factor  which  the  HVOSM  had  to  be  able  to 
dupl icate. 

The  test  vehicle,  a  Vega,  was  well  documented  so  that  the  innut 
parameters  were  assured  to  match  the  test  vehicle  (2).  Figure  7  lists 
the  card  images  of  the  input  for  the  modified  HVOSM  simulation.  The 
variables  associated  with  the  input  numbers  in  Figure  7  can  be  found  in 
the  HVOSM  User's  Manual  (3).  The  accelerometer  readings  from  the  right 
side  of  the  vehicle  were  used  as  direct  approximations  of  the  horizontal 
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Figure  7.  Input  for  Crash  Test  Comparisons. 


31 


forces  generated  in  the  collision.  This  was  because  the  signpost  im¬ 
pacted  15  inches  to  the  right  of  center,  so  that  accelerometer  was 
closest  to  the  lateral  point  of  contact.  Then,  since  the  post  had  bent 
around  the  hood  of  the  car  and  imparted  large  vertical  forces  to  the 
front  of  the  car,  the  vertical  forces  were  successively  approximated 
until  the  HVOSM  results  matched  the  crash  test  results.  The  final 
applied  vertical  impulse  appears  in  Figure  8.  The  point  of  application 
for  the  horizontal  and  vertical  impulses  was  15  inches  to  the  right  of 
center  on  the  front  bumper,  just  as  in  the  crash  test.  Also,  the  accel¬ 
erometer  placements  on  the  simulated  HVOSM  vehicle  matched  those  on  the 
actual  crash  test  vehicle. 

The  results  of  the  simulation  are  shown  in  Table  1  and  Figure  9  . 
Using  all  of  the  available  recorded  data  to  check  the  accuracy  of  the 
modified  HVOSM  program  simulation,  the  program  results  matched  ex¬ 
tremely  well  with  the  crash  test  #11  results.  As  shown  in  Table  1,  the 
maximum  error  was  below  10%  (except  for  isolated  accelerometer  points), 
and  the  accelerometer  readings,  both  taken  from  the  right  accelerometer 
of  the  vehicle  shown  in  Figure  9,  correlate  very  closely.  The  simu¬ 
lated  vehicle  rolled,  as  did  the  test  vehicle,  and  at  almost  the  same 
time.  Based  on  relatively  little  data,  the  modified  HVOSM  was  able  to 
duplicate  almost  exactly  the  results  of  crash  test  #11.  Just  to  prove 
that  this  was  not  an  isolated  success,  a  second  crash  test  exercise  was 
chosen  as  a  final  challenge  to  the  acceptance  of  the  modified  HVOSM. 

Crash  Test  Number  9 

Crash  test  #9  (7),  as  performed  by  TTI,  was  similar  to  crash  test 
#11  in  many  ways.  Crash  test  #9  used  the  same  Vega,  used  the  same 
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accelerometer  placement,  and  had  the  same  type  of  recorded  data  as 
crash  test  #11.  There  were  some  challenging  differences,  however. 
First  of  all,  the  impact  occurred  15  inches  to  the  left  of  center,  so 
the  left  accelerometer  readouts  were  used  to  simulate  the  horizontal 
impulse  forces  with  time.  Second  of  all,  the  signpost  was  markedly 
different  and  much  stiffer  than  AASHTO  standards.  While  the  signpost 
in  crash  test  #11  was  not  AASHTO-approved,  the  signpost  in  crash  test 
#9  was  far  in  excess  of  AASHTO  standards,  and  as  such  represented  an 
extreme  test.  Finally,  the  automobile  nearly  rolled,  but  it  did  not 
experience  rollover  even  though  the  signpost  wrapped  around  the  vehi¬ 
cle  front  as  in  crash  test  #11.  Therefore,  the  challenge  was  to  take 
the  car  to  the  brink  of  rollover  but  not  over  it,  while  duplicating 
the  rest  of  the  test  results  as  closely  as  possible.  Using  succes¬ 
sive  approximations,  as  in  the  crash  test  #11  duplication,  the  final 
vertical  sign  impulse  input  developed  for  the  exercise  is  shown  in 
Figure  10. 

As  seen  in  Table  2  and  Figure  11,  the  modified  HVOSM  came  even 
closer  to  matching  the  actual  crash  test  results  than  in  the  previous 
exercise.  The  errors  were  more  than  halved  from  what  they  were  in 
the  crash  test  #11  effort,  and  the  accelerometer  readouts  in  Figure 
11  once  again  shadowed  those  of  the  actual  crash  test.  The  vehicle 
lost  ground  contact  on  three  of  the  four  wheels,  but  did  not  roll 
over,  just  as  had  been  planned.  Another  result,  not  shown  in  the 
table  or  on  a  graph,  was  that  it  took  only  one-third  of  the  number  of 
iterations  used  in  the  exercise  for  crash  test  #11  to  duplicate  crash 
test  #9.  In  other  words,  experience  made  the  validation  exercise 
easier  to  complete. 


Vertical 
Crash  Tes 
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Another  interesting  sidelight  can  be  observed  by  comparing  the 
signpost  vertical  impulses  for  the  two  crash  tests,  shown  in  Figure  8 
(p.  32)  and  Figure  10  (p.  36).  Two  different  signs,  both  of  which  did 
not  break  away  but  which  bent  around  the  hood  and  caused  large  vertical 
impulses,  almost  matched  vertical  impulse  curves  in  separate  develop¬ 
ment  of  each  curve.  The  vertical  impulse  of  the  crash  test  #11  curve 
was  547  pound-seconds,  while  the  vertical  impulse  of  the  crash  test  #9 
curve  was  561  pound-seconds,  a  spread  of  only  1.3".  This  coincidence 
was  used  later  to  define  the  vertical  impulse  in  the  parameter  study. 

Conclusion 

The  modified  HVOSM  had  to  prove  its  validity  by  comparison  with 
both  theory  and  actual  signpost  impact  test  results.  Completing  each 
challenge,  the  modified  HVOSM  was  shown  to  be  reliable  and  easy  to 
use.  The  error  in  many  cases  was  negligible,  and  the  duplications 
were  accurate  for  even  the  most  complex  tests.  Using  relatively 
sparse  amounts  of  crash  test  data,  the  modified  HVOSM  was  able  to 
duplicate  two  different  signpost  collisions  almost  exactly,  matching 
all  measured  output  from  the  real  crash  tests.  With  the  completion 
of  this  phase  of  the  thesis,  the  parameter  study  of  signpost  impacts 
could  be  realistically  and  authoritatively  examined. 


CHAPTER  IV 


PARAMETER  STUOY  APPLICATION 

Introduction 

After  proving  the  modified  HVOSM  to  be  an  accurate  simulator  of 
signpost  collisions,  the  crux  of  the  thesis  could  be  properly  ap¬ 
proached.  AASHTO  standards,  as  mentioned  at  the  start,  were  developed 
before  the  trend  toward  subcompact  cars  became  evident.  Since  no  para¬ 
meter  studies  comparing  the  effects  of  different  sized  cars  in  colli¬ 
sions  with  a  signpost  meeting  the  upper  limit  of  AASHTO  standards  had 
been  performed,  the  HVOSM  was  modified  for  the  purpose  of  just  such  a 
parameter  study. 

Using  the  same  horizontal  and  vertical  impulses  on  each  size  car, 
all  approaching  at  the  same  speed  before  the  collision,  the  modified 
HVOSM  was  used  to  compare  the  effects  of  these  impulses  on  the  differ¬ 
ent  automobiles.  Two  different  signpost  types  were  used:  the  break¬ 
away  model,  which  consisted  of  only  horizontal  impulses;  and  the  base¬ 
bending  model,  which  consisted  of  both  horizontal  and  vertical  impulses. 
The  automobile  parameters  were  input  for  a  large  car,  a  medium  car,  a 
compact  car,  and  a  subcompact  car.  The  credibility  of  the  parameter 
study  can  be  better  shown  by  examining  the  preparation  of  the  sign  and 
automobile  models  upon  which  the  study  results  were  based. 

Signpost  Models 


As  mentioned  above,  two  signpost  models  were  used  in  the  study. 
The  models  were  developed  using  AASHTO  standards  (10)  and  the  results 
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of  the  validation  tests.  The  AASHTO  standards,  in  brief,  limit  the 
change  in  momentum  to  1100  pound-seconds  for  a  2250  pound  automobile 
striking  a  support  at  speeds  between  20  and  60  miles  per  hour.  These 
specifications  apply  for  both  breakaway  signposts,  which  are  designed 
to  break  away  at  ground  level,  and  base-bending  signposts,  which  bend 
at  the  base  as  the  vehicle  rides  down  the  length  of  the  signpost. 

These  two  signpost  models  required  impulses  to  fit  the  AASHTO  specifi¬ 
cations. 

The  breakaway  model  signpost  was  composed  of  only  horizontal  im¬ 
pulses.  The  horizontal  impulse,  shown  in  Figure  12,  was  developed 
using  the  left  accelerometer  readout  from  crash  test  #11  (7),  and  the 
total  impulse  is  1095  pound-seconds,  just  less  than  the  upper  limit  set 
by  AASHTO.  Examining  Figure  12  in  detail,  the  total  impulse  up  to 
.185  seconds  is  1095  pound-seconds,  with  the  subsequent  impulses 
negating  each  other  arithmetically.  That  .185  second  impulse  value  is 
important  because  it  is  the  time  point  used  to  calculate  the  change  in 
momentum  in  each  program  simulation.  This  signpost  impulse  was  used 
for  the  breakaway  model  impact  simulation  input  with  each  vehicle  size, 
from  large  to  subcompact. 

The  base-bending  model  signpost  was  composed  of  both  horizontal 
and  vertical  impulses.  The  horizontal  impulse  used  was  the  same  as 
that  used  in  the  breakaway  model.  This  choice  was  made  assuming  that 
the  same  signpost  used  in  the  breakaway  model  was  driven  into  the 
ground  rather  than  mounted  to  break  away  at  the  base.  The  vertical 
impulse,  produced  as  the  signpost  wrapped  around  the  front  of  the 
vehicle,  was  taken  as  a  combination  of  the  vertical  impulses  found  in 
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the  duplication  of  crash  tests  #11  and  #9.  As  previously  noted,  the 
impulses  of  those  crash  tests  were  547  pound-seconds  and  561  pound- 
seconds  for  tests  #11  and  #9,  respectively.  Since  these  two  different 
posts  applied  almost  identical  vertical  impulses  in  base-bending  sign¬ 
post  impacts,  the  vertical  impulse  for  the  base-bending  model  was  based 
on  these  impulses.  The  resultant  vertical  impulse,  shown  in  Figure  13, 
measured  567  pound-seconds.  The  small  error  is  due  to  the  fact  that 
the  time  points  had  to  be  those  used  for  the  horizontal  impulse  input, 
since  the  horizontal  impulse  was  most  critical.  The  input,  after  ad¬ 
justing  to  the  horizontal  impulse  time  steps,  contained  some  roundoff 
errors  in  the  vertical  impulse  which  accounted  for  the  small  impulse 
increase. 

In  short,  the  breakaway  model  and  the  base-bending  model  were 
based  on  AASHTO  standards  and  actual  crash  test  results.  The  break¬ 
away  model  consisted  of  a  horizontal  impulse  only,  while  the  base¬ 
bending  model  consisted  of  both  horizontal  and  vertical  impulses  to 
simulate  the  effects  of  a  signpost  wrapping  around  the  front  of  the 
vehicle  as  the  car  rides  down  the  signpost.  Both  models  were  developed 
for  impact  by  the  four  car  sizes  to  examine  how  the  same  signposts 
could  affect  different  sizes  of  automobiles  in  a  collision.  With  the 
completion  of  the  signpost  models,  the  only  required  input  development 
was  the  car  data  for  the  operation  of  the  HVOSM  simulations. 

Vehicle  Data 

Four  vehicle  sizes  were  used  in  this  parameter  study:  large, 
medium,  compact,  and  subcompact.  Each  of  these  had  to  be  defined  as 
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accurately  as  possible  for  the  program  input  variables.  In  some  cases, 
the  car  data  were  incomplete,  so  that  other  data  which  would  least 
affect  the  results  were  substituted.  For  example,  the  medium  and  sub¬ 
compact  cars  lacked  some  minor  program  values,  so  the  values  for  the 
large  and  compact  cars,  respectively,  were  used  for  those  gaps  in  the 
data.  It  was  felt  that  using  the  data  from  the  next  largest  vehicle 
than  the  one  lacking  data  would  at  least  serve  to  prevent  exaggeration 
of  program  results.  To  illustrate,  using  large  car  suspension  data  for 
the  medium  car  would  moderate  the  results  to  a  degree,  while  using  the 
compact  vehicle's  suspension  data  might  magnify  the  impact  effects  on 
the  medium  car.  Also,  each  vehicle  had  the  impulse  positioned  on  the 
front  bumper  15  inches  to  the  right  of  center.  Since  the  front  bumper 
is  approximately  at  the  same  height  as  the  center  of  gravity,  no  dis¬ 
placement  along  the  vehicle  z-axis  was  input  for  any  car,  but  the  dis¬ 
tance  to  the  front  bumper  measured  along  the  vehicle  x-axis  from  the 
center  of  gravity  varied  with  the  car  size. 

Several  other  parameters,  none  associated  with  a  particular  car, 
were  input  in  all  cases.  The  road  surface  coefficient  of  friction  was 
set  at  .5  for  all  tests,  even  though  many  tire  manufacturers  claim  a 
friction  factor  of  .8.  This  decision  was  made  based  on  the  fact  that 
signposts  are  placed  off  the  road  shoulder  and  not  on  a  paved  surface. 
All  vehicles  approached  impact  at  60  miles  per  hour,  corresponding  to 
the  upper  limit  set  in  the  AASHTO  standards.  All  tests  were  performed 
using  the  free-steering  mode  of  the  modified  HVOSM  operation  to  study 
the  free-wheel  stability  of  the  vehicle  in  post-impact  motion.  Finally, 
all  program  simulations  were  run  for  3  seconds,  which  was  a  long  enough 
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time  to  observe  the  important  accelerations  and  vehicle  motions. 

Large  car.  One  of  the  two  well -documented  vehicles,  the  large  car 
was  represented  by  a  1963  Ford  sedan  (2).  All  vehicle  data  necessary 
for  the  modified  HVOSM  were  available,  and  the  input  as  loaded  for  the 
program  appears  in  Figure  14.  The  data  set  ICARD  values  appear  as  the 
last  number  in  every  line.  As  can  be  seen  by  examining  the  ICARD=29 
series,  the  collision  input  data  is  for  the  breakaway  model  signpost  as 
distinguished  by  the  lack  of  z -forces  in  the  data.  As  explained  in  the 
crash  test  duplication  efforts,  the  variables  associated  with  each 
number  can  be  found  using  the  HVOSM  User's  Manual  (3).  An  explanation 
of  each  variable  as  presented  in  the  input  would  be  tedious  and  detract 
from  the  main  points  of  the  parameter  study. 

Medium  car.  While  the  large  car  input  was  handled  easily,  the 
data  for  the  medium  car  posed  some  problems.  No  complete  listing  of 
the  required  variables  was  available  for  the  medium  car.  After  a 
search,  the  best  source  was  found  in  a  compilation  of  vehicle  data  for 
parameter  studies  (1).  This  listed  many  of  the  required  variable 
values  for  a  general  medium  car,  composed  of  the  average  values  for 
many  medium  cars  and  thus  not  associated  with  a  make  or  model  year. 
While  incomplete,  the  data  were  the  best  which  could  be  found,  and  all 
gaps  were  filled  using  the  large  car  data.  The  important  and  most 
influential  variables  for  the  study,  however,  were  supplied  in  the 
medium  car  data;  and  the  large  car  data  were  required  mainly  to  fill 
out  the  required  input.  As  shown  in  Figure  15,  the  medium  car  input 
was  therefore  more  nearly  correct  than  the  inclusion  of  large  car  data 
might  make  it  seem.  To  emphasize,  the  large  car  data  were  used  for 


IN»UT  ttCAD  BV  CAL SVA 


46 


o 

iTt 

43 


<rNi?-0‘C‘^(?o(rao^ 

—  xrv,  t\J  <\  <v  r'r\  <v<v  <m  <V 


O 


??Of  Pffff  P  ’  r>  -  -  3  •■»  3  >  ■■»  .'(??* 
P  (<♦*'  ( -r*C>r  '•  *  *  Mf\J^r 

a 

j 


o  •  o  oo 

•  IT  O  Iflo 

x  l*.*  oo 

i  <» «  •  f- 


o 

o 

£ 


o 


{ 

i 


•  <V 

(V.  Ifl 

o 

OOO  •  •  I/’  o 

o 

9  ON  • 

r. 

©, 

M  I  X  •  -O 

• 

o 

•  •  • 

o  c « 
ooo*  c 

o 

C  ♦ 

O  c  o 

**  <o  9  tr  c 

o 

o 

F*  1  Oalfl  • 

o 

# 

• 

•  <M 

o 

OO* 

O  Og  X  X  (T  O 

1 

• 

4>  *OON  • 

• 

o 

«-CO  •©  • 

« 

•  <r  rf*  r  &  i**  *r  ,p  it  f>  ui 

N 

r  -  »  *-  o 

i 

o 

c 

• 

o<\»  o  o 

M 

ni 

o  *  •  «  o  • 

• 

1 

a 

o  •©pk>o  • 

be  500  )0  <300  3 

• 

0303^  O 

*■ 

•oflutrtTttf 

* 

o 

<0 

o 

a 

<t  0  • 

/I 

o 

•41  41N-  3 

o 

3 

O  ‘ION  %o  3 

o 

oooooooooa© 

• 

V  O  •  •  •  •  o 

• 

<\* 

— 

n  4)  xm«  a  x 

o 

• 

xOOOOJCOOOJ 

OcOOCCoCC  ©c  OCCOCQOOC 

i 

i 


*r  «•  *o  kP it  w> *r  »'  f  -r  n in  >r -r t~  jiio •* 


3  0  300  3  0'3CHJ-03'3oCli'3J3 
tab  cr:  s  r  r  rcr.:  sjfr  cH 


C0©033>0  =  V  -»  3  J  C  3  000  0  00 


930i33000d  J  ^OOoO.-^-iQOi- 


CWCoO'-.  OC  -  CiOC'.r  a  c  c.  -*>  n  - 
oooooooo  •  o  LoO:  .,',03 


•  •  •  • 

«”4kl  •••  •  • 

o  **  m  o  •  •  M  t  -£  ■  •  a  •  tfl  •  ■ 

y  <.A  c  r  •*  •  «  r-  p;  4  tcro^n 

C~«»''''f'*xf'iO*~^''>Nfr,K  03000 

—  O  i  rt  -  %  • 

I  I  t  t  .-<  1^  .*  I  •  I  I  I  *  I  W  '«00 


rj  J>  •  •  n 

On  ■••«>> 

•  *A*n  jo 
-  J  .flx..- 


oo  • 
•  *m 
OO  I 


•  <i  oj  o  «  r>  <r 
oo**ojmr»^iri4SP‘ic 
o  oooOooooo 
. . . 


x  *-  m  <r-  -n  '*  *<.  r  <t  —  ,r  —  r  r  o  *  *  x  ►» 

>Joo«','*'fl3KNn:>oo«Mw  • 


Figure  14.  Input  for  Large  Car. 


1NPUI  MIAU  bY 


47 


W 

C 


—  —  <ii  <\<  .n  ;v  ru  py  cii  ru  <M  r«  o* 


>3>0-993*0' 

PvJ^rX<YfY«\i<'JM<M'MvM-MiVj'W<MfNjryfv»!')^' 


O  «3 
•  03 

-'no 
i  n* 


UO 

wo 

0*Y 

•  w 


o 

3W  • 

oTl 


o 

•  wo 

*n  • 

in  •  o 


•  i 

•  m  a 

a  -w 
<v  on  © 
O  — WW*^  • 
•  «3  • 

O- I -- • 


/)**•--  9 

**>•<3  o  OPY  •  • 

ON  *0  3  *9  •« 

N  —  ®  •  •  M  O  I 


W*—  09 

\  C  CO  •  •  o  •  • 

34  • flN  30  «W 

•  o  — 


♦  • 

•  ntn  o 

o  o  •  <>*-  — o«© 
•  a  <o  *  •  •  •  o  • 
o^wnwwo-o 


•  oo  • 

9  •  •O  • 
-OO  I  - 


333 

30 030000000 


. . 

•*nwww%-nnwwwww 

oooaJcocoooaooc 

• 

O 

9 

O  «•••••• 

O  •  *030  3  0  0  0 

oo  oo  owooooo  o 

o  •  «wo»no  —  n<j«cY 

-33M4-«lfl  «  *  *  W 


NO  0090000000 
W  •»••••••••• 

•  O  oo  00000009 


•  •  S3  •  •  •  —  »K  — 

ooNon«w*N»fYO 

9  .<?--< 

—  O  I 


00003000000 30003000 


w  o  w  w  w  w  w  w  w  w  W  V»  W  W  w  VI  w  w  o 


1}Q4}49oaOXO(0(O^Q  O  ®  <X  343 

N*».f**.  f»M»  h  - 


rfffrfr: 


300993030303 
3  011(50  33350C3T33  30  30 
OO  M  N  W  W  n  3  3N  0*1  •••«••• 
O0Yi0rtrt^r\N\.liN330  30  30 


003303900003003903  3 
•  ••«••••<•(••••••«• 

003300  3  30  3333039033 

•  •  • 

cww  •  •  •  •  * 

%  —  W  •  •P4<fC»*0«v,'«  • 

vl3«0  •  WN*Y*T3CO«iO  •  O 

—  — no  «  -N^s^NMOCt 9000 

— •  • 

t  j  I  »  i  i  l  *  3  I  (  I  NO  «03 


W(*g  •‘JJ  <J  •  W<V  ONW^WW  — 

*n  •  •>»>  oo  •  -  oo-NAPifdo^S? 

•  — O  •  »w  I  O  *0003031533 

0*3--  ®Ot  1  ^3  •  I  «  «  •  l  •  l  i  « 


3>  t  -N 

>  3  O  n  n  *  W^)N  ■»*  -5  3  ^  3  3  —  ^  <*J  , 

3  —  —  —  —  -J  ‘g  JOi'JC 


» 


Figure  15.  Input  for  Medium  Car. 
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minor  variables  which  were  not  covered  in  the  referenced  work. 

An  examination  of  Figure  15  also  reveals  that  the  impulse  input 
card  series  is  that  for  the  base-bending  signpost  model,  highlighted  by 
the  inclusion  of  vertical  as  well  as  horizontal  forces.  Also,  com¬ 
paring  Figure  14  (p.  46)  and  Figure  15  will  show  that  the  bumper  posi¬ 
tion  has  moved  closer  to  the  vehicle  center  of  gravity  from  the  large 
to  the  medium  car,  as  expected. 

Compact  car.  Another  well -documented  vehicle,  the  compact  car  in¬ 
put  data  was  easily  compiled  (2).  Shown  in  Figure  16,  the  data  is  that 
for  a  1971  Chevrolet  Vega  sport  coupe,  the  same  basic  vehicle  as  that 
used  for  the  crash  tests  covered  previously.  All  variable  values  had 
been  determined,  so  input  presented  no  problem.  Figure  16  also  pre¬ 
sents  the  base-bending  signpost  impulse  data,  and  a  quick  comparison 
with  Figure  15  will  support  the  fact  that  the  impulse  data  is  identi¬ 
cal  in  each  signpost  model,  regardless  of  vehicle  size. 

Subcompact  car.  Presenting  a  final  problem,  the  subcompact  car 
data  was  the  hardest  to  find.  Eventually,  a  source  containing  the  most 
important  vehicle  parameters  was  found  (9).  As  with  the  medium  car, 
the  rest  of  the  required  input  data  had  to  be  obtained  from  the  next 
largest  vehicle.  Although  more  Vega  data  were  required  for  the  sub¬ 
compact  (a  1976  Honda  Civic)  than  had  been  required  for  the  medium  car, 
the  crucial  and  most  influential  data  were  found  for  the  Honda,  so  the 
Vega  values  were  merely  supplemental  and  had  little,  if  any,  effect  on 
the  final  results.  Figure  17  is  a  listing  of  the  subcompact  input  for 
the  HVOSM,  also  showing  the  same  breakaway  model  impulse  data  as  pre¬ 
sented  in  Figure  14  (p.  46)  for  the  large  car.  As  before,  the  bumper 
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Figure  16.  Input  for  Compact  Car. 
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Figure  17.  Input  for  Subcompact  Car. 
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continued  to  move  closer  along  the  vehicle  x-axis  to  the  vehicle  center 
of  gravity  as  the  car  size  decreased,  but  the  15  inches  right  of  center 
offset  of  the  force  application  was  maintained. 

Conclusion 

Approaching  the  parameter  study  with  the  necessity  of  credible 
results  in  mind,  the  input  had  to  be  reliable  for  the  output  to  present 
a  realistic  approximation  of  actual  signpost  impacts.  The  signpost  im¬ 
pulse  data  sets  were  developed,  one  for  a  breakaway  model  signpost,  and 
the  other  for  a  base-bending  signpost  model.  The  input  for  the  vehicle 
data  was  developed  from  the  best  available  data  for  each  of  the  car 
sizes:  large,  medium,  compact,  and  subcompact.  With  the  preparation 
of  the  input  data,  each  of  the  four  vehicle  sizes  were  prepared  for 
simulated  impact  with  each  of  the  signpost  models. 

At  this  point,  the  modified  HVOSM  was  put  into  operation  to  com¬ 
pute  the  results  of  each  of  the  eight  separate  vehicle  cases.  Analysis 
and  presentation  of  those  results  was  the  final  step  of  the  thesis. 
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CHAPTER  V 

PARAMETER  STUDY  RESULTS 

Introduction 

Modifying  the  HVOSM,  validating  the  changed  program,  and  assem¬ 
bling  data  necessary  to  evaluate  AASHTO  signpost  standards  culminated 
in  the  results  shown  in  this  chapter:  the  reactions  of  different  sized 
automobiles  in  collisions  with  each  of  two  different  signpost  configu¬ 
rations.  Points  of  interest  are  not  only  whether  the  effect  of  the 
same  signpost  is  markedly  different  on  decreasing  sizes  of  vehicles, 
but  also  how  different  methods  of  ground  connection  for  the  signpost 
(breakaway  versus  base  bending)  affect  the  same  vehicle  in  a  collision. 

Parameter  study  results  will  be  delineated  in  three  major  cate¬ 
gories:  the  longitudinal  acceleration  of  the  sprung  mass  center  of 
gravity,  the  vertical  acceleration  of  the  sprung  mass  center  of  gravity, 
and  the  vehicle  travel  path  after  impact.  Center  of  gravity  accelera¬ 
tions  are  a  measure  of  the  accelerations  felt  by  the  occupants  of  the 
vehicle,  while  the  vehicle  travel  path  is  an  indication  of  post-impact 
stability  of  the  vehicle. 

Before  the  separate  results  are  presented,  a  quick  view  of  the 
overall  results  will  help  to  keep  the  major  trends  in  focus.  First  of 
all,  as  expected,  the  changes  in  momentum  were  all  in  accordance  with 
AASHTO  specifications  with  one  exception:  the  subcompact  vehicle,  when 
impacted  with  the  base-bending  signpost  model,  measured  1385  pound- 
seconds  of  momentum  change,  26%  above  AASHTO  standards.  In  this 
instance,  the  subcompact  rolled  over,  even  though  the  same  impulse  did 
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not  affect  the  larger  three  cars,  in  terms  of  change  in  momentum  or  in 
terms  of  vehicle  motion,  nearly  as  severely.  As  anticipated  at  the  out¬ 
set  of  the  thesis,  the  subcompact  experienced  extreme  accelerations  and 
one  rollover  in  tests  with  the  two  signpost  models. 

Longitudinal  Acceleration 

The  longitudinal  acceleration  of  the  sprung  mass  center  of  gravity 
was  measured  as  an  indication  of  the  front-to-rear  accelerations  expe¬ 
rienced  by  the  vehicle  passengers.  The  accelerations  were  measured 
only  for  the  duration  of  the  impact  (up  to  .25  seconds),  since  they 
were  negligible  after  that  period  of  time.  Generally,  the  maximum 
deceleration  for  each  vehicle  size  was  only  slightly  less  for  the  base¬ 
bending  model  than  for  the  breakaway  model.  Also,  the  maximum  decel¬ 
eration,  as  expected,  increased  for  each  decreasing  vehicle  size,  al¬ 
most  doubling  between  the  compact  and  subcompact  vehicle  sign  impact 
simulations. 

Large  car.  The  large  car  accelerations  are  shown  in  Figure  18. 

The  maximum  deceleration  is  only  3.45  G's,  experienced  in  the  breakaway 
model  of  signpost  impact.  The  values  for  the  two  different  signpost 
models  are  in  almost  exact  agreement  for  the  first  .16  seconds.  After 
that,  however,  the  acceleration  values  are  small  and  the  differences 
between  the  acceleration  curves  would  not  have  an  appreciable  affect  on 
the  passengers,  apart  from  the  slightly  greater  decelerations  expe¬ 
rienced  in  the  base-bending  model. 

Medium  car.  The  medium  car  accelerations  are  shown  in  Figure  19. 
Graphed  to  the  same  scale  as  the  large  car  accelerations  shown  in 


Breakaway  Model 


Longitudinal  Acceleration 
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Figure  18,  the  maximum  negative  acceleration  has  increased  to  5  G's, 
again  for  the  breakaway  model  of  sign  impact.  As  before,  the  accelera¬ 
tions  are  in  close  agreement  for  the  two  sign  impact  models,  separating 
only  after  .13  seconds  into  the  collision.  As  with  the  large  car,  the 
base-bending  model  of  signpost  collision  ends  up  causing  greater  decel¬ 
erations  than  the  breakaway  model  at  the  later  stages  of  the  collision. 
While  the  shape  of  the  acceleration  curve  is  similar  to  that  for  the 
large  car,  the  magnitudes  are  greater  overall  for  the  medium  car,  and 
the  larger  accelerations  and  rapid  reversals  would  be  more  noticeable 
to  the  vehicle  passengers. 

Compact  car.  Presented  in  Figure  20,  the  compact  car  accelerations 
are  greater  than  those  for  the  medium  car  at  all  points,  much  like  a 
larger  version  of  the  same  graph.  Increasing  to  5.49  G's,  the  maximum 
negative  acceleration  is  experienced  in  simulated  impact  with  the  break¬ 
away  model  signpost.  As  before,  the  breakaway  and  base-bending  models 
match  accelerations  for  .13  seconds,  after  which  the  base-bending  model 
is  associated  with  greater  deceleration  than  the  breakaway  model.  The 
accelerations,  while  greater  in  magnitude  than  those  for  the  medium  car, 
are  not  relatively  that  much  larger  than  those  for  the  medium  car,  so 
the  effects  on  the  vehicle  would  be  only  slightly  greater. 

Subcompact  car.  The  differences  between  the  accelerations  of  the 
subcompact  car  and  the  three  larger  vehicles  are  of  such  magnitude  that 
the  scale  of  the  graph  had  to  be  increased.  Shown  in  Figure  21,  the 
acceleration  scale  is  double  the  scale  that  was  used  for  the  three 
larger  vehicles.  While  it  appears  to  be  the  same  as  that  for  the  com¬ 
pact  car,  the  graph  of  the  accelerations  of  the  subcompact  car  is 
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actually  almost  twice  the  magnitude.  Another  significant  change  is 
found  in  the  difference  between  the  accelerations  experienced  under  the 
two  sign  models  in  Figure  21.  Starting  at  .09  seconds  after  impact, 
the  base-bending  signpost  model  causes  much  greater  decelerations  by 
the  vehicle  than  the  breakaway  model  does.  A  maximum  negative  acceler¬ 
ation  of  9.59  G's,  as  experienced  in  the  breakaway  model,  is  almost 
double  the  maximum  of  the  compact  car  and  almost  triple  that  of  the 
large  car.  While  all  accelerations  are  for  short  durations,  the  passen¬ 
gers  in  the  subcompact  car  would  be  in  much  greater  danger  from  decel¬ 
eration  effects  than  they  would  be  in  any  of  the  three  larger  vehicles. 
As  previously  mentioned,  rollover  occurs  in  the  base-bending  model 
simulation  for  the  subcompact,  but  the  high  accelerations  experienced 
in  the  breakaway  model,  even  without  rollover,  are  hazardous. 

In  brief,  the  front- to-rear  accelerations  of  the  vehicles  increase 
in  magnitude  as  the  car  size  decreases.  While  all  accelerations  are  of 
short  duration,  they  are  so  markedly  greater  for  the  subcompact  car  as 
to  render  its  collision  effects  on  the  vehicle  passengers  particularly 
undesirable  relative  to  the  accelerations  experienced  in  any  of  the 
three  larger  vehicles.  Also,  the  base-bending  model  accelerations 
agree  closely  with  the  breakaway  model  accelerations,  differing  after  a 
period  of  time  in  the  direction  of  larger  decelerations  of  the  vehicle. 
As  the  car  size  decreases,  the  deviancy  between  the  sign  model  accel¬ 
erations  grows  in  magnitude  and  the  deviancy  begins  a  shorter  time 
after  impact. 
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Vertical  Acceleration 

At  the  sprung  mass  center  of  gravity,  the  vertical  accelerations 
were  compiled  as  a  measure  of  the  up-and-down  accelerations  felt  by  the 
passengers  of  each  vehicle  during  signpost  impact.  People  have  lower 
tolerances  to  head-to-toe  accelerations  than  they  have  to  chest-to-back 
accelerations,  which  are  longitudinal  accelerations  in  this  study. 
Therefore,  the  vertical  accelerations  in  these  cases  have  as  much  im¬ 
portance  in  measuring  vehicle  safety  as  the  longitudinal  accelerations. 
These  vertical  accelerations  were  graphed  for  a  period  of  .25  seconds, 
which  was  long  enough  to  cover  all  the  extreme  acceleration  values  asso¬ 
ciated  with  the  impacts.  In  general,  while  the  breakaway  signpost  only 
produced  vertical  accelerations  worthy  of  notice  in  the  subcompact  car 
test,  the  base-bending  model  signpost  produced  high  vertical  accelera¬ 
tions,  especially  in  the  smaller  vehicles. 

Large  car.  Shown  in  Figure  22,  the  graph  of  the  accelerations  ex¬ 
perienced  by  the  sprung  mass  center  of  gravity  does  not  resemble  the 
impulse  input  as  the  longitudinal  accelerations  resembled  the  hori¬ 
zontal  impulse  input.  There  is  no  graph  of  the  accelerations  produced 
in  the  breakaway  model  because  the  maximum  acceleration  of  .033  G’s 
would  not  show  up  on  the  graph.  An  interesting  reaction  is  produced  by 
the  initial  positive  acceleration  of  the  vehicle  center  of  gravity, 
followed  by  high  negative  accelerations.  Since  the  positive  direction 
of  acceleration  is  downward  for  the  HVOSM  vehicle  (Appendix  II),  that 
means  that  the  passengers  would  feel  the  seats  drop  away,  then  slam 
upwards  and  drop  away  again  as  the  negative  acceleration  peak  at  .215 
seconds  occurs.  While  not  critical  in  this  case  due  to  the  constructive 


62 


effects  of  seatbelts  and  the  almost  instantaneous  existence  of  the 
accelerations,  these  positive  and  negative  accelerations  grow  in  magni¬ 
tude  while  acting  in  the  same  short  time  span  as  the  car  size  decreases. 
As  will  be  shown,  the  greater  magnitudes  and  large  reversals  point  to 
increasingly  hazardous  conditions,  even  with  seatbelts  and  shoulder 
harnesses. 

Medium  car.  Once  again,  the  breakaway  model  fails  to  produce  any 
significant  vertical  accelerations,  with  the  maximum  acceleration  only 
reaching  .09  G's.  However,  in  the  base-bending  model  signpost  impact 
simulation,  the  maximum  negative  acceleration  increases  over  60%  from 
that  of  the  large  car  impact,  as  shown  in  Figure  23.  As  before,  the 
initial  downward  acceleration  is  followed  by  rapid  negative  accelera¬ 
tion  (upward),  which  then  quickly  drops  off.  These  rapid  reversals  are 
the  real  danger  since  the  other  vehicle  motions  are  not  very  extreme 
(as  will  be  shown).  Driver  control  becomes  more  difficult  as  these 
bounding  motions  become  more  violent. 

Compact  car.  The  sprung  mass  center  of  gravity  vertical  accelera¬ 
tions  as  recorded  for  the  base-bending  model  sign  impact  are  shown  in 
Figure  24.  At  first  glance  the  effects  seem  to  be  reduced,  but  the 
further  increase  in  acceleration  magnitudes  caused  the  scale  of  accel¬ 
erations  to  be  doubled.  So,  the  accelerations  are  really  shown  to  half 
size  of  what  they  would  be  on  the  scale  used  for  the  medium  and  large 
cars.  The  usual  downward  acceleration  followed  by  the  large  upward 
acceleration  creates  the  same  dangerous  situation  for  the  automobile 
occupants,  and  the  maximum  upward  acceleration  increased  10%  over  the 
maximum  for  the  medium  car.  As  with  the  two  larger  vehicles,  there  are 
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23.  Sprung  Mass  C.G.  Vertical  Acceleration 
for  Medium  Car  in  Base  Bending  Model. 
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not  any  significant  vertical  accelerations  produced  in  the  breakaway 
model  signpost  impact. 

Subcompact  car.  Shown  on  the  same  scale  as  those  for  the  compact 
car  (making  the  accelerations  appear  only  half  as  large  as  they  would 
be  on  the  scale  used  for  the  medium  and  large  cars),  the  accelerations 
of  the  vehicle  under  both  signpost  models  appear  in  Figure  25.  The 
breakaway  model  produces  larger  vertical  accelerations  for  the  subcom¬ 
pact  car  than  for  any  of  the  other  car  sizes,  but  even  the  accelerations 
experienced  in  this  case  are  negligible  and  do  not  present  a  hazard, 
especially  relative  to  any  of  the  accelerations  experienced  in  the 
base-bending  model  of  signpost  impact.  For  an  idea  of  the  magnitude  of 
the  maximum  negative  acceleration  for  the  subcompact  in  the  base¬ 
bending  model,  the  value  is  3.5  times  the  maximum  experienced  by  the 
large  car  and  twice  as  large  as  that  for  the  compact  car.  The  base¬ 
bending  model,  it  must  be  remembered,  causes  rollover  of  the  subcompact, 
but  the  high  vertical  accelerations  would  probably  render  the  vehicle 
occupants  incapable  of  noticing  the  vehicle  rollover.  The  vertical 
accelerations  caused  by  the  base-bending  model  on  the  subcompact  are 
even  greater  than  the  maximum  longitudinal  acceleration  (for  which  the 
body  has  a  higher  tolerance). 

As  was  the  case  with  the  longitudinal  accelerations,  the  vertical 
accelerations  of  the  vehicle  increase  in  magnitude  as  the  vehicle  size 
decreases.  Under  the  base-bending  model  impact,  the  maximum  vertical 
acceleration  is  greater  than  the  maximum  longitudinal  acceleration. 
Unfortunately,  the  vertical  accelerations  affect  the  body  more  than  the 
longitudinal  accelerations.  The  breakaway  model  signposts  produce  the 
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Figure  25.  Sprung  Mass  C.G.  Vertical  Acceleration 
for  Subcompact  Car. 
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greatest  longitudinal  accelerations,  but  the  model  only  produces 
noticeable  vertical  accelerations  on  the  subcompact  vehicle. 

Travel  Path 

The  lateral  motion  of  the  vehicle  under  a  free  steering  mode  gives 
an  indication  of  the  controllability  of  the  vehicle  after  sign  impact. 
The  following  graphs  show  both  the  path  of  the  vehicle  and  an  arrow 
at  intervals  to  give  an  idea  of  the  vehicle  orientation  as  it  travels 
after  the  signpost  impact. 

In  presenting  the  vehicle  paths,  full  use  of  each  graph  is  made  to 
present  the  information  in  an  uncluttered  manner.  The  side  effect  of 
this  decision,  however,  is  a  great  variance  in  scale  of  the  lateral 
displacement  of  the  vehicle.  Therefore,  the  analysis  below  is  made 
with  an  effort  to  keep  the  relative  paths  of  each  vehicle  in  perspec¬ 
tive.  Mainly,  though,  the  problem  is  that  the  graphs  cannot  be  com¬ 
pared  without  adjusting  for  scale  when  viewing  them. 

Vehicle  stability  after  collision  is  measured  by  the  possible  con¬ 
trollability  of  the  vehicle  after  signpost  impact.  For  example,  a  car 
which,  in  the  free  steering  mode,  corrects  itself  smoothly  to  its  pre¬ 
impact  orientation  is  very  stable;  on  the  other  hand,  a  wildly  spinning 
car  could  not  be  contolled  by  any  driver  input,  and  thus  the  vehicle 
reaction  is  classified  as  unstable.  Of  course,  any  car  which  rolls 
over  does  not  require  a  path  analysis;  the  hazard  is  obvious. 

Large  car.  Performing  as  expected,  the  large  car  is  very  stable 
in  either  sign  model  impact  simulation.  Shown  in  Figure  26,  the  break¬ 
away  model  sign  impact  results  in  a  10-foot  displacement  to  the  right 
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by  the  large  car.  Travelling  smoothly,  the  car  fishtails  very  slightly 
to  the  left  at  sign  impact,  correcting  its  own  course  during  travel 
with  subsequent  steer  to  the  left.  On  the  same  scale,  the  base-bending 
model  sign  impact  vehicle  motion  is  presented  in  Figure  27.  Small 
steering  angles  result  in  only  a  6-foot  displacement  to  the  right. 

Both  vehicle  paths  could  be  easily  controlled  by  any  driver  input,  and 
the  lateral  travel  distance  of  the  vehicle  relative  to  the  "forward" 
travel  of  the  vehicle  is  very  small. 

Medium  car.  Presented  on  two  different  lateral  scales,  the  travel 
paths  of  the  medium  car  show  good  stability.  After  impact  with  the 
breakaway  signpost,  the  medium  car  steers  wide  right  and  then  starts  to 
overcorrect  to  the  left  as  illustrated  in  Figure  28.  While  travelling 
24  feet  to  the  right,  the  vehicle  can  be  easily  controlled  before  the 
overcorrection  to  the  left  by  the  free  steering  vehicle  takes  place. 

The  base-bending  signpost  impact  travel  path,  presented  in 
Figure  29,  is  even  more  stable  and  has  less  lateral  travel.  The  lat¬ 
eral  travel  of  only  9  feet,  compared  to  24  feet  for  the  breakaway  sign¬ 
post,  places  the  travel  path  of  the  vehicle  for  the  base-bending  model 
impact  very  close  to  the  original  course  compared  to  the  breakaway 
model  impact.  However,  a  growing  factor  in  the  driver  controllability 
is  the  effect  of  the  longitudinal  and  vertical  accelerations  previously 
mentioned  on  the  driver's  ability  to  steer  the  car.  While  neither  of 
the  medium  car  impact  models  would  incapacitate  the  driver,  the  fol¬ 
lowing  two  car  sizes  demonstrate  more  hazardous  outcomes,  especially 
in  the  base-bending  model  impacts. 

Compact  car.  As  just  stated,  the  high  accelerations  experienced 
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in  the  vehicle  collisions  with  the  signpost  models  could  have  affected 
the  driver's  ability  to  steer  the  vehicle.  But,  just  as  the  vehicle 
stability  is  getting  more  critical  due  to  increasing  driver  hazard  from 
accelerations,  the  vehicle  stability  decreases  with  decreasing  vehicle 
size.  As  shown  in  Figure  30,  the  compact  car  impacts  with  the  break¬ 
away  sign  and  fishtails  gradually  left  until  the  vehicle  is  oriented 
at  right  angles  to  its  original  course.  A  lateral  displacement  of  50 
feet  is  unsafe  because  other  obstacles,  such  as  fences,  ditches,  and 
embankments,  exist  along  highways;  and  many  of  these  obstacles  are 
within  50  feet  of  the  road  shoulder.  While  the  essentially  longitu¬ 
dinal  accelerations  generated  in  the  breakaway  model  impact  would  not 
incapacitate  the  compact  car  driver,  the  car  is  seen  as  marginally 
stable  because  of  the  excessive  lateral  travel  and  the  need  for  control. 

The  base-bending  model  impact,  on  the  other  hand,  produces  a  rela¬ 
tively  small  lateral  displacement  of  17  feet  as  seen  in  Figure  31.  The 
vehicle  corrects  its  own  path,  which  is  important  because  the  vertical 
and  longitudinal  accelerations  would  almost  certainly  limit  the  driver’s 
ability  to  control  the  path  of  the  vehicle.  As  mentioned  at  the  start, 
the  scale  differences  between  the  travel  path  graphs  along  the  lateral 
axis  for  the  compact  car  are  very  great,  and  the  lateral  displacement 
of  the  compact  car  impacting  the  base-bending  model  signpost  would  be 
very  small  on  the  graph  of  the  travel  path  for  the  breakaway  signpost 
impact. 

Subcompact  car.  Complete  loss  of  stability  is  demonstrated  in 
each  signpost  model  collision  simulation  using  the  subcompact  car.  The 
breakaway  model  causes  vehicle  spinout,  while  the  base-bending  model 
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figure  30.  Path  and  Orientation  of  Compact  Car 
in  Breakaway  Model. 
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impact  results  in  vehicle  rollover.  No  matter  which  signpost  model  is 
used,  the  vehicle  path,  when  combined  with  the  vehicle  accelerations, 
prints  only  to  extreme  hazard.  Presented  in  Figure  32,  the  path  and 
orientation  of  the  subcompact  car  subsequent  to  impact  with  the  break¬ 
away  model  signpost  reveals  the  360+°  spin  and  the  relatively  small 
lateral  displacement  of  only  8  feet.  That  is  the  smallest  breakaway 
model  sign  impact  lateral  displacement  of  the  four  vehicle  sizes,  and 
the  main  reason  for  the  small  lateral  displacement  is  the  rapid  spin  of 
the  vehicle.  No  vehicle  control  is  possible  for  the  subcompact,  but, 
after  the  large  longitudinal  accelerations  experienced  by  the  driver 
in  the  sign  impact,  there  is  not  a  high  expectation  of  driver  input 
after  impact. 

The  limited  results  of  the  base-bending  model  impact  are  shown  in 
Figure  33.  The  rollover  is  considered  completed  after  a  45°  roll  angle 
is  established.  The  rapid  spinning  of  the  vehicle  prior  to  rollover  is 
shown,  along  with  the  fractional  lateral  displacement.  Effectively,  in 
comparison  with  all  other  vehicle  paths,  the  travel  path  prior  to  roll¬ 
over  for  the  subcompact  in  the  base-bending  model  impact  was  quite 
similar  to  that  of  the  large  car  in  the  base-bending  model  impact. 

In  general,  the  stability  of  the  vehicle  after  signpost  impact  de¬ 
creases  with  decreasing  vehicle  size,  illustrated  in  this  study  by  the 
total  loss  of  control  of  the  subcompact  car.  The  base-bending  model 
signpost  impact  tends  to  result  in  more  stable  motion  and  less  lateral 
displacement  than  the  breakaway  sign  impact.  The  increased  stability 
of  the  base-bending  model  impact  is  offset  by  the  higher  longitudinal 
and  vertical  accelerations  endured  by  the  vehicle  occupants. 
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Summary 

A  compilation  of  the  essential  parameters  of  the  study  are  gathered 
in  Table  3.  The  different  car  sizes  and  their  results  in  each  of  the 
two  signpost  model  simulated  impacts  are  presented  for  easy  comparison 
and  for  review  of  the  trends  covered  in  this  chapter. 

First  of  all,  the  accelerations  as  measured  at  the  sprung  mass 
center  of  gravity,  and  thus  approximately  the  same  as  those  experienced 
by  the  vehicle  occupants,  increase  as  the  car  size  decreases.  The 
breakaway  model  signpost  is  associated  with  slightly  larger  longitudinal 
accelerations  and  insignificant  vertical  accelerations  compared  to  those 
recorded  in  the  base-bending  model  impacts.  In  each  case,  however,  the 
maximum  accelerations  suffered  by  the  subcompact  are  approximately 
double  those  of  the  compact  car  and  triple  those  of  the  large  car. 

Secondly,  the  vehicle  grows  more  unstable  as  the  vehicle  size  de¬ 
creases;  this  is  compounded  by  the  increased  incapacity  of  the  driver 
due  to  the  increasing  accelerations  with  decreasing  vehicle  size.  The 
subcompact  is  totally  unstable,  experiencing  spinout  in  the  breakaway 
signpost  impact  and  rollover  in  the  base-bending  model  signpost  impact. 
The  same  signpost  that  produces  small  pitch  and  yaw  angles  (shown  in 
Table  3)  in  the  compact  vehicle  is  devastating  to  the  subcompact  in  the 
base-bending  model  impact  simulations.  While  the  vehicle  occupants 
were  safer  in  the  breakaway  model  impact  of  the  subcompact,  the  sur¬ 
prising  fact  is  this:  The  difference  between  installing  a  breakaway 
base  support  and  simply  driving  the  same  signpost  into  the  ground  (base 
bending)  in  a  collision  by  a  subcompact  vehicle  is  the  difference  be¬ 
tween  a  survivable  spinout  and  a  possibly  fatal  vehicle  rollover. 


TABLE  3.  Vehicle  Sign  Impact  Results 
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CHAPTER  VI 
CONCLUSION 

As  times  change,  standards  must  be  kept  up  to  date  with  technolo¬ 
gical  breakthroughs  and  changing  conditions  which  affect  those  stan¬ 
dards.  The  recent  trend  in  smaller,  lighter  cars  brought  about  by  the 
gas  shortage  and  the  use  of  plastics  technology  to  decrease  vehicle 
weight  warranted  the  re-examination  of  signpost  standards,  especially 
in  light  of  crash  tests  which  indicated  a  hazard  for  light  vehicles  in 
signpost  impacts  (7).  While  the  standards  for  signposts  had  been  esta¬ 
blished  in  consideration  of  compact  cars  when  compact  cars  were  a  small 
minority  of  road  traffic,  the  growing  percentages  of  compact  and  sub¬ 
compact  cars  make  roadside  design  for  the  smaller  sizes  a  necessary 
concern. 

To  properly  study  the  effects  of  a  sign  deemed  safe  by  AASHTO 
standards  on  a  range  of  car  sizes,  a  dependable  and  valid  computer 
simulation  had  to  be  used  in  order  to  obtain  repeatable  results.  The 
best  model,  the  HVOSM,  was  modified  and  validated  to  accept  idealized 
sign  impact  data  as  generalized  impulses.  The  three-dimensional  capa¬ 
bilities  of  the  modified  HVOSM  provided  an  excellent  means  to  study 
the  post-impact  motions  and  their  effects  on  the  vehicle.  Four  vehicle 
sizes  were  used  to  cover  the  general  vehicle  categories:  large,  medium, 
compact,  and  subcompact.  Two  signpost  models  were  created  to  repre¬ 
sent  the  general  installation  methods:  breakaway  at  the  base,  and 
base  bending  due  to  placement  of  the  signpost  into  the  ground  as  its 
own  support. 

The  input  for  each  vehicle  and  signpost  model  was  made  with 
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deliberate  care  to  approximate  real  impact  conditions.  The  impulses 
were  generated  to  test  the  extremes  of  the  AASHTO  standards;  that  is, 
a  60  mile  per  hour  impact  with  a  signpost  impacting  an  impulse  of  1100 
pound-seconds.  The  breakaway  signpost  was  idealized  as  causing  only 
horizontal  impulses,  while  the  base-bending  model  contained  the  same 
horizontal  impulse  plus  a  vertical  impulse  modeled  after  actual  crash 
test  results. 

The  results  of  the  simulation  confirm  the  expected  results:  The 
signposts  which  are  safe  for  the  larger  cars  have  hazardous  effects  on 
the  subcompact  car.  The  breakaway  sign  impact  results  in  spinout  and 
excessively  high  accelerations  for  the  subcompact,  and  the  base-bending 
model  impact  causes  even  greater  accelerations  and  eventual  vehicle 
rollover.  In  each  case,  the  base  bending  model  of  installation  pro¬ 
duces  more  hazardous  effects  for  the  vehicle  occupants  than  the  same 
signpost  installed  for  breakaway  action  in  impacts. 

While  these  results  definitely  point  to  the  need  for  re-examination 
of  the  AASHTO  standards,  they  also  point  to  an  extension  of  these  para¬ 
meter  studies.  This  parameter  study  was  performed  assuming  ideal 
roadside  conditions:  flat,  level  terrain  with  a  coefficient  of  fric¬ 
tion  between  the  tires  and  the  surface  of  .5.  Realistic  roadside  condi¬ 
tions  would  tend  to  magnify  the  hazards  of  impact  with  sideslope, 
shoulder  roughness,  and  possible  roadside  ditches.  The  proposed  new 
areas  of  investigation  could  be  easily  approached  using  the  modified 
HVOSM,  a  proven  simulator  of  signpost  impact  effects.  Regardless  of 
further  studies,  this  thesis  reveals  that  the  AASHTO  standards  expose 
drivers  of  subcompact  vehicles  to  dangerous  sign  impact  effects  even 
under  the  best  of  conditions. 
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APPENDIX  II 

A  DESCRIPTION  OF  THE  HIGHWAY-VEHICLE-OBJECT 
SIMULATION  MODEL 

The  Highway-Vehicle-Object  Simulation  Model  (HVOSM)  is  a  mathema¬ 
tical  model  which  can  be  utilized  to  investigate  various  problems  asso¬ 
ciated  with  the  roadway  environment.  It  is  used  to  study  the  dynamic 
behavior  of  an  automobile  under  different  roadway  and  roadside  condi¬ 
tions.  These  conditions  include  such  problems  with  the  roadway  envi¬ 
ronment  as  highway  traffic  barrier  collisions,  rapid  lane  change  maneu¬ 
vers,  handling  response  on  horizontal  curves,  drainage  ditch  cross 
sections,  and  now  sign  impacts. 

The  HVOSM  was  developed  by  Cornell  Aeronautical  Laboratory  (CAL) 
(5,6)  and  later  modified  for  specific  program  studies  by  the  Texas 
Transportation  Institute  (TTI)  (11).  A  conceptual  idealization  of  the 
model  is  shown  in  Figure  34.  The  model  is  idealized  as  four  rigid 
masses,  which  include:  (a)  the  sprung  mass  (Ms)  of  the  body  supported 
by  the  springs,  (b)  the  unsprung  masses  (Mi  and  M2)  of  the  left  and 
right  independent  system  of  the  front  wheels,  and  (c)  the  unsprung  mass 
(M3)  representing  the  rear  axle  assembly. 

The  eleven  degrees  of  freedom  of  the  model  include  translation  of 
the  automobile  ’ n  three  directions  measured  relative  to  the  fixed 
global  coordinate  system,  rotation  about  the  three  local  axes  of  the 
automobile,  independent  displacement  of  each  front  wheel  suspension 
system,  suspension  displacement  and  rotation  of  the  rear  axle  assembly, 
and  steer  of  the  front  wheels.  A  further  explanation  of  this  model  for 
in-depth  analysis  can  be  found  by  examining  the  references  quoted 
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earlier  in  this  section,  and  also  by  examining  Reference  8. 
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APPENDIX  III 

CORRECTION  OF  ORIGINAL  HVOSM  FOR  IBM  360  COMPUTER 
AND  CURRENT  FORTRAN  H  EXTENDED  COMPILER 


Introduction 

The  HVOSM  program  deck,  which  was  the  1976  version  stored  in  the 
files  at  Texas  A&M  University  as  supplied  by  CALSVA  and  modified  by 
TTI,  would  not  run  without  program  error  statements  generated  by  the 
Fortran  H  extended  compiler.  The  correction  of  these  errors  was  very 
time  consuming  and  frustrating,  so  this  documentation  will  list  the 
changes  made  in  order  to  get  the  version  to  run  prior  to  modification. 
The  errors  were  of  three  basic  types:  subroutine  name  errors,  EQUI¬ 
VALENCE  statement  errors,  and  BLOCK  DATA  errors. 

Subroutine  Name  Errors 

When  the  HVOSM  was  written,  the  CLEAR  subroutine,  which  is  built 
into  the  computer  system  and  not  listed  in  the  card  deck,  was  used  in 
the  main  program  and  the  MATRIX  subroutine  to  clear  addresses  in  the 
common  blocks.  The  assumed  operation  of  subroutine  CLEAR  is  that  the 
parameter  list  contains  the  first  and  last  addresses  in  the  common 
block,  and  the  computer  clears  all  addresses  between,  and  including, 
the  two  addresses  given  in  the  parameter  list.  The  program  would  not 
run  on  the  Texas  A&M  computer  system  because  this  subroutine,  while  it 
is  listed  in  the  computer  system,  is  not  the  same  one  that  the  devel¬ 
opers  of  the  HVOSM  used. 

On  the  Texas  A&M  IBM  360,  the  subroutine  ERASE  is  provided  which 
uses  the  same  parameter  list  and  performs  the  same  function  that  the 
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CLEAR  subroutine  did  in  the  original  program  on  whatever  computer  sys¬ 
tem  CALSVA  used.  All  that  had  to  be  done  to  correct  this  error  was  to 
replace  the  CLEAR  subroutine  call  statements  in  MAIN  and  MATRIX  with 
ERASE  subroutine  call  statements  —  a  process  of  merely  replacing  the 
word  CLEAR  with  the  word  ERASE. 

EQUIVALENCE  Statement  Errors 

This  was  by  far  the  largest  source  of  errors  in  the  HVOSM  program. 
When  the  program  was  originally  written,  the  compiler  would  perform 
different  functions  than  the  present  compiler  in  use  at  Texas  A&M.  The 
explanation  of  the  reason  for  changes  and  the  list  of  necessary  changes 
is  very  simple,  but  the  actual  program  changes  were  very  numerous  and 
time  consuming.  This  explanation  is  intended  to  save  a  lot  of  time, 
and  the  action  to  clear  up  these  errors  results  in  a  much  neater  and 
more  compact  program. 

Using  the  older  compiler  systems,  an  equivalence  could  be  made 
between  a  matrix  and  a  certain  series  of  variables  quite  easily.  For 
example,  the  matrix  A(4)  and  the  variables  Al,  A2,  A3,  and  A4  could  use 
the  following  EQUIVALENCE  statement: 

EQUIVALENCE(A.Al) 

This  statement  would  automatically  set  A(l)  equal  to  Al,  A(2)  equal  to 
A2,  and  so  on  until  the  limits  of  the  matrix  were  met.  This  operation 
was  possible  because  the  EQUIVALENCE  statement  took  the  first  variable 
of  a  numerical  sequence  and  set  the  equivalence  between  the  first 
matrix  location  and  the  variable  listed,  and  then  automatically  set  the 
equivalence  for  the  same  number  of  following  addresses  in  the  coronon 
block  as  there  were  in  the  matrix.  In  another  example,  two  equally 
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sized  matrices  could  be  set  equal  to  each  other  with  a  single  EQUIVA¬ 
LENCE  statement.  If  the  matrices  were  named  ONE  and  TWO,  the  computer 
would  equate  the  values  in  similar  locations  of  each  matrix  so  that 
they  were  identical  by  the  use  of  the  following  command: 

EQU I VALENCE ( ONE , TWO ) 

The  computer  compilers  at  Texas  A&M,  however,  do  not  allow  these 
two  operations.  Each  EQUIVALENCE  statement  must  have  a  one-to-one 
correspondence  between  variables  listed  in  the  statement,  which  re¬ 
quires  the  filling  of  each  separate  location  of  matrices.  So,  in  the 
past  it  was  easier  to  link  several  different  subroutines  and  their 
individual  variables  using  EQUIVALENCE  statements  and  matrices.  Now, 
however,  it  is  easier  and  also  neater  to  scrap  all  duplicate  variables 
so  that  one  program  variable  is  used  throughout  the  program. 

The  following  list  of  EQUIVALENCE  statements  which  had  to  be  re¬ 
paired  are  accompanied  by  an  explanation  of  the  program  surgery  used 
to  repair  the  errors.  All  of  the  listed  EQUIVALENCE  statements  were 
deleted  from  the  program  after  they  were  replaced  by  the  new  EQUIVA¬ 
LENCE  statements  or  they  were  rendered  useless  by  the  elimination  of 
one  of  the  variables  in  each  set.  All  changes  were  made  to  the  pro¬ 
gram  to  affect  the  least  number  of  equations,  and  all  "series"  varia¬ 
bles  like  Al,  A2,  A3,  and  A4  were  replaced  by  the  simpler  matrix  form. 

EQUIVALENCE(YCIP.YCIP) 

The  change  was:  YC1P  was  replaced  by  YCIP(l),  and  YC2P  was  re¬ 
placed  by  YCIP(2).  The  common  block  INPT1  was  changed  by  replacing 
"YC1P,  YC2P,"  with  "YCIP(2) ,"  and  the  DIMENSION  statement  for  YCIP  was 
then  deleted.  Then,  all  the  YCIP  and  YC2P  locations  were  replaced  by 
YCIP(l)  and  YCIP(2),  respectively. 


EQUIVALENCE^, X1P),(YP, YIP),  (ZP, ZIP), (PHII.PHIl), 
(PSII  ,PSI1)  ,(UI  ,U1)  ,(VI  ,V1) ,  (WI  ,V/1) 


90 


This  statement  was  corrected  by  the  lengthy  procedure  of  replacing 
X1P  through  X4P  with  XP(1)  through  XP(4)  throughout  the  whole  program. 
All  of  the  XP,  YP,  ZP,  PHI I ,  PSII,  U,  V,  and  W  matrices  are  one¬ 
dimensional  and  four  addresses  long.  So,  in  like  manner,  YIP  to  Y4P, 

ZIP  to  Z4P,  PH II  to  PHI4,  PSII  to  PSI4,  U1  to  U4,  VI  to  V4,  and  Wl  to 
W4  were  replaced  by  YP(4),  ZP{4),  PHII(4),  PSII(4),  UI(4),  VI(4),  and 
WI(4),  respectively,  with  the  variable  numbers  matching  the  matrix 
address  locations  in  each  case. 

The  common  blocks  were  also  changed  to  reflect  these  improvements. 
DIMV  common  block  was  modified  so  that  XP(4),  YP(4),  ZP(4),  PHI  1(4), 
and  PSII (4)  replaced  the  previously  listings  of  the  variables  X1P  to 
X4P,  YIP  to  Y4P,  and  so  on.  All  of  the  equations  in  which  any  of  the 
phased  out  variable  names  were  listed  were  changed  to  reflect  the  new 
matrix  names  and  locations.  The  common  block  ADTNL  was  modified  so 
that  UI(4),  VI(4),  and  WI(4)  replaced  the  previous  listings  of  U1  to  U4, 
VI  to  V4,  and  Wl  to  W4.  Those  equations  had  to  be  changed  for  the  new 
matrix  names  and  locations,  also. 

Special  caution  had  to  be  used  in  subroutine  DRIVER,  where  the 
variable  PSII  was  used  in  the  parameter  list.  In  the  DRIVER  subroutine, 
the  variable  names  for  PSII  were  changed  to  PSII,  but  were  not  sub¬ 
scripted  with  a  matrix  location  because  the  parameter  list  would  not 
accept  a  matrix  location,  and  also  because  it  was  a  local  subroutine 
variable  and  not  used  in  the  normal  program  sense  in  that  subroutine. 

After  all  of  this  was  performed,  the  DIMENSION  statement  which  had 
originally  declared  the  matrices  was  then  deleted  as  redundant,  since 
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the  common  blocks  DIMV  and  ADTNL  now  handled  that  function. 

EQUIVALENCE(INDXPT,I1) 

This  statement  was  corrected  by  replacing  II  through  14  with 
INDXPT(l)  through  INDXPT(4),  respectively.  This  involved  extensive 
modification  of  equations  in  the  AREA,  IDOUT,  and  SFORCE  subroutines. 

In  particular,  the  IDOUT  subroutine  presented  problems  in  statement 
numbers  IDOT  187  and  IDOT  262.  In  those  lines,  there  are  implied  DO 
loops  in  the  WRITE  statements.  The  index  value,  II,  of  each  of  the 
DO  loops  is  bounded  by  the  initial  value  II  and  the  test  value  14. 

These  needed  to  be  replaced  by  INDXPT(l)  and  IN0XPT(4 ) ,  but  the  array 
locations  could  not  be  used  in  the  implied  DO  loop  as  the  initial  and 
test  values.  Instead,  the  addition  of  the  following  two  statements 
immediately  preceding  the  affected  WRITE  statements  resolved  the 
problem: 

IPIGZ=INDXPT(1) 

I DAWGZ= I NDXPT ( 4 ) 

These  two  integer- valued  variables  were  then  inserted  into  each  WRITE 
statement  implied  DO  loop  as  the  initial  and  test  values,  respectively. 

EQUIVALENCE(XGP,THG) , (YGP, PHIG) 

This  pair  of  statements  had  to  be  corrected  by  replacement,  since 
all  of  the  listed  variables  are  21  x  21  size  matrices  and  a  listing  of 
the  equivalence  of  each  location  would  have  been  too  lengthy.  After 
checking  all  of  the  listings  of  the  individual  variables,  the  logical 
replacement  choice  was  to  replace  THG  and  PHIG  by  XGP  and  YGP,  res¬ 
pectively.  This  was  decided  because  there  were  many  new  equations  with 
XGP  and  YGP,  and  the  fact  that,  although  THG  and  PHIG  were  listed  in 
all  the  INPT  comnon  blocks,  all  the  cotnnon  block  replacements  could  be 
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made  with  one  command  due  to  the  unique  listing  of  the  variables,  sepa¬ 
rated  by  commas,  which  is  found  only  in  the  common  blocks  and  other 
deceleration  statements.  Therefore,  using  the  WYLBUR  terminal,  one 
CHANGE  command  eliminated  most  of  the  locations  of  THG  and  PHIG. 

The  only  other  location  where  THG  and  PHIG  were  listed  was  in  the 
BLOCK  DATA  section  of  the  program,  where  simple  substitution  in  the 
Hollerith  string  was  made.  This  same  line  was  included  as  one  of  the 
BLOCK  DATA  errors  covered  later  in  this  appendix. 

In  conclusion,  the  variable  matrices  THG  and  PHIG  were  replaced  by 
XGP  and  YGP,  which  involved  a  simple  substitution  in  the  INPT  common 
block  and  one  line  in  the  BLOCK  DATA  section  of  the  HVOSM. 

EQUIVALENCE ( YPNO , P ) 

EQU I VALENCE (YN0,Q) 

YPNO  and  YNO  were  two  one-dimensional  arrays,  30  addresses  lone. 
These  variables  were  all  confined  to  the  PINT1  subroutine,  and  YPNO  and 
YNO  were  the  least  used  pair  of  variables  out  of  each  pair  of  possible 
variable  choices.  Since  there  were  so  few  references  to  YPNO  and  YNO 
in  the  subroutine,  it  was  a  simple  matter  to  replace  them  with  P  and  Q, 
respectively.  Both  EQUIVALENCE  statements  were  deleted,  and  YPNO  and 
YNO  had  to  be  deleted  from  a  DIMENSION  statement  and  a  DOUBLE  PRECISION 
statement  in  the  PINT1  subroutine  to  complete  this  correction. 

EQU I VALENCE ( HCAH , HCAH 1 ) , ( HCBH , HCBH1 ) , ( HCGH , HCGH 1 ) 

This  EQUIVALENCE  statement  was  corrected  in  a  different  way  from 
all  of  the  previous  ones.  There  were  considerable  numbers  of  each 
variable  in  the  program,  so  that  no  clear  choice  presented  itself  as  in 
the  other  EQUIVALENCE  statements.  Also,  this  EQUIVALENCE  statement, 
along  with  the  arrays  HCAH,  HCBH,  and  HCGH  only  occurred  in  the  VGORNT 
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subroutine.  Since  this  statement  only  caused  an  error  in  one  subrou¬ 
tine,  and  since  each  array  was  only  a  one-dimensional,  four  location 
array,  the  solution  used  in  this  case  was  to  list  each  location  in  each 
array  to  an  equivalent  variable.  Therefore,  the  EQUIVALENCE  statement 
was  enlarged  so  that  each  of  the  four  locations  in  the  arrays  HCAH, 
HCBH,  and  HCGH  were  linked  to  the  variables  HCAH1  through  HCAH4,  HCBH1 
through  HCBH4,  and  HCGHI  through  HCGH4,  respectively.  As  previously 
stated,  this  was  a  unique  case  because  of  the  great  number  of  listings 
of  each  variable  in  the  program  and  the  fact  that  this  substitution 
only  had  to  be  performed  in  the  VGORNT  subroutine. 

BLOCK  DATA  Errors 

There  were  three  simple  errors  in  the  BLOCK  DATA  section  of  the 
program,  all  of  which  were  easily  corrected  and  which  are  mentioned  as 
a  timesaver  in  case  the  corrections  need  to  be  made  again. 

First  of  all,  in  line  BDAT  7,  there  was  a  closure  of  the  data  set 
listing  before  the  array  was  filled.  The  BAX  data  set  was  defined  as 
30  locations  long,  but  the  string  ran  out  after  22  places  were  filled. 
Therefore,  at  the  end  of  the  string,  eight  spaces  had  to  be  filled  as 
follows: 

»8*‘  7 

This  will  fill  up  the  remaining  eight  spaces  in  the  array  and  then 
eliminate  the  error  message. 

The  second  error  was  almost  identical  to  the  first.  The  BEX  array 
was  also  defined  as  30  locations  long,  but  the  data  listed  only  21 
spaces  before  it  ended.  Therefore,  the  following  was  added  after  the 
last  parenthetical  listing: 


This  added  the  required  nine  spaces  to  the  string  to  complete  the  array 
and  eliminate  the  error  message. 

The  final  BLOCK  DATA  error  message  was  due  to  Hollerith  string 
errors.  The  original  listing,  in  part,  was  shown  in  subroutine  IDOUT 
as: 

DATA  IDDATA/12HZPG,THG,PHIG/, — 

The  mistake  here  is  that  three  variables  have  been  defined  in  one 
Hollerith  string,  with  the  commas  included.  As  mentioned  in  the 
EQUIVALENCE  error  portion,  THG  and  PHIG  were  replaced  by  XPG  and  YPG, 
respectively,  so  the  correct  listing  ends  up  as: 

DATA  I DDATA/3HZPG , 3HXPG , 3H YPG/ , — 
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APPENDIX  IV 

PROGRAM  MODIFICATION  LISTING 

This  additional  section  lists  the  computer  printouts  of  the  changes 
made  to  the  HVOSM  in  order  to  add  the  sign  impact  option,  along  with 
the  regular  program  statements  surrounding  each  modified  area  so  that 
the  precise  location  of  the  changes  in  the  HVOSM  can  be  found. 

Figure  35  is  a  listing  of  the  changes  made  to  the  main  program. 
While  it  does  not  include  the  placement  of  COLL,  SUBB,  or  NUMT  in  the 
common  blocks,  that  was  explained  sufficiently  well  that  it  was  not 
required.  As  the  last  variable  in  the  BARIER  common  block,  the  ERASE 
statement  had  to  be  changed  also  to  clear  COLL  and  SUBB  along  with  the 
rest  of  the  conmon  block  by  inserting  SUBB(7,50)  in  the  place  of  the 
second  argument  in  the  parameter  list  for  that  program  statement. 

Figure  36,  covering  three  pages,  illustrates  the  changes  made  to 
the  INPUT  subroutine,  while  Figure  37  illustrates  the  two  additions 
made  to  the  SFORCE  subroutine  at  the  beginning  and  end  of  the  sub¬ 
routine.  Finally,  Figure  38  lists  the  formats  and  variable  lists  used 
to  identify  and  arrange  the  printout  of  the  SUBB  matrix  in  the  IDOUT 
subroutine. 
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COMPUTED  TO  IT*IT+1  AFTER  RETURNIN6  FROM  SVAPLT,  gHICH  ALLOWS  OPEN  ENDERAIN  113 


CONCERNING  DATA  SETS  AND  SEQUENTIAL  DATA  SETS  PER  JOB  FOR  PRIVATE  PAK  NAIN  11A 

1T=8  RAIN  115 

CET  NPTS  FOR  NUflBER  OF  WRITES  =  NURBER  POINTS/CATEGORY (RAIN  TO  PLTPJRAIN  116 

t;******************************************************************************* 

C***  The  CORRON/BARIER/  BLOCK  HAS  BEEN  RUBIFIED  FOR  THE  SIBN  IRPACT  ****** 

C*«*  VERSION  UF  THIS  PROGRAR  BY  THE  ADDITION  OF  THE  'COLL'  TR16GER  ****** 

t*«*  OPT 1UN  VALUE  AND  THE  'SU88!7,50)'  RATRIX  OF  IRPACT  FORCES  AND  ****** 

c*«*  their  locations  uith  tire,  me  currand  below  this  clears  the  ****** 

c***  VALUE  OF  'COLL-  FOR  LATER  IN  THE  PROGRAR.  THE  CORRON/IRPT1/  ****** 

L’***  BLOCK  HAS  BEEN  ROOIFIED  BY  THE  ADDITION  OF  THE  'HURT'  VARIABLE  ****** 

C***  WHICH  WILL  BE  USED  TO  INDICATE  THE  RURBER  OF  TIRE  POINTS  USED  ****** 

C***  TO  READ  IN  THE  IRPACT  FORCES  AND  LOCATIONS.  IT  IS  ALSO  CLEARED  ****** 

C***  TO  A  VALUE  OF  ZERO  BELOW  AT  THE  START  OF  THE  PROGRAR.  ****** 

t*************** ***************** t*********4*»***********«**********»*********** 

NUflT=0 

NPTS=1  RAIN  117 

BEX(30)=0.  RAIN  118 

1  CALL  ERASE ! PHI O.PQRRIN)  RAIN  11* 

CALL  ERASE! YC1P! 1 ) , PSIFOO )  RAIN  120 


Figure  35 


Modification  of  MAIN  Program 


L**« 

L'*** 

L*** 

C*** 

C*** 

L*** 


DATA  SAVE/ 'SAVE'/ 
DATA  FINl/'flNI'/ 
DATA  STN  /'STAN'/ 
BATA  TNUNE  /'NONE'/ 
NSAVf»0 


INPT  1M 
I  HPT  115 
INPT  116 
INPT  117 
- INPT  TT8 


THIS  IS  PANT  Of  THE  SI8N  1NPACT  NUDIFICATION  FUR  THE  CURRENT  ****** 

HVOSN  US1N6  AN  INPACT  NATRIX  OF  THE  FORCES  AND  LOCATIONS  UITH  TINE  ****** 

THE  NEXT  LINE  CLEARS  A  VARIABLE  USED  TO  LOAD  THE  HATRIX  FROH  ****** 

****** 


THE  INPUT  DATA 


JJ*0 

IAPFRO )  a  0 
1APFR<2)  *  0 
JUT*0 
ALPHA=0 . 0 


INPT  119 
INPT  120 
INPT  121 
INPT  122 


2  READ  1001,  (CARDIH(K),K*1,20),(DUNU),i*l,9) , I CARD  INPT  133 

1001  FORNAT <  20 A*, I1,9F8.0,I8)  INPT  136 

UR I T£ <6,101 1 )  (CARD1N(K),K=1 ,20)  INPT  135 

1011  FORNATt  1X.20A4)  INPT  136 

IFUCARD.0E.y999)  GOTO  60  INPT  137 


C******************************************************************************* 

L*«*  THE  'UO  TO'  OPTIONS  LISTED  BELOW  HAS  BEEN  NODIFIED  BY  THE  ADDITION  ****** 

C***  OF  THE  LABEL  53  FOR  WHEN  ICARD  EOUALS  29,  TO  READ  IN  THE  INPACT  ****** 

l***  FORCES  AND  LOCATIONS  WITH  TINE  AS  NODIFIED  BY  JOHN  BETZ  ****** 

(.'*•******•***•*******••••**••*•****•**•**•********•**•**•**•*•**•**••*•*•***«*»* 
60  TO  (3, A, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16,19, 20, 21, 22, 23, 24, 26, 32, INPT  138 
1  34,35,36, 37,38,39,53) , ICARD  INPT  139 

WRITE(6,2001 )  (CARDIN(K),K*1 ,20)  INPT  140 


Figure  36.  Modification  of  INPUT  Subroutine 
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in  *  BUIU3)  MPT  223 

AAA  *  DUflM  >  INPT  22A 

HAAX  *  DUfl<5)  I APT  225 

HN1N  =  D0N<6)  INPT  226 

BET  *  DUN<7)  INPT  227 

C******************************************************************************* 

C**«  THIS  IS  PAKT  Of  THE  VEHICLE  iflPACT  BUDIF1CAHON  PERFQRflED  ON  THE  ****** 

CM*  PR06RAB.  THE  VARIABLE  'COLL'  HAS  FEEfT  ADDETTAS' A  SIGNAL  USED  ****** 

CM*  LATER  IN  THE  SfORCE  SUBROUTINE  TO  PREVENT  THE  BARRIER  OPTION  ****** 

C«*«  FROH  C0HIN6  INTO  PLAY  WITHOUT  AVUIDINS  THE  REST  OF  THE  PROGRAN.  ****** 

C***  THIS  ALSO  ALLOWS  THE  BARRIER  OPTION  TO  BE  KEPT  AND  THEREFORE  THE  ****** 

C«**  PROORAN  IS  AVAILABLE  FOR  A  VARIETY  OF  PROBLENS.  THE  OTHER  ****** 

C*»*  VARIABLE,  'NOHT',  IS  THE  READ-IN  VALUE  OF  THE  NUNBER  OF  ****** 

C*«*  DISCRETE  TINE  STEPS  WHICH  WILL  BE  READ  INTO  THE  PROORAN  TO  ****** 

C***  REPRESENT  THE  IflPACT  FORCES  AND  LOCATIONS  FOR  THE  SINULATION.  ****** 

CM*  THE  VALOE  OF  'NOflT'  HOST  BE  BETWEEN  1  AND  50,  AND  IT  SHOULD  ****** 

CM*  INCLUDE  a  ZERO  TIflE  VALUE,  THE  FINAL  TINE  VALUE  WHEN  ALL  THE  ****** 

C**«  FORCES  ARE  ZERO,  AND  THE  TIflE  POINT  IN  THE  FUTURE  WHICH  WILL  KEEP  ****** 

C***  THE  FORCES  ZE.H’OED  OUT  FOR  THE  PROSRAfl  AS  IT  KEEPS  RUNNIN6.  THE  ****** 

cm*  the  last  tiae  point  should  be  a  long  period  between  the  next-to-  ****** 

C*M  LAST  SO  THE  PROORAN  CAN  RUN  THE  VEHICLE  REACTION.  ****** 

;;  ******************************  ********  *****•***«••  *****************  ******  «***«* 


C0LL=DUfl(8) 
NUNTaDON( P  > 
GO  TO  2 
5  XNS  *  DUflUJ 
XNUF  =  DUN<2) 
xnuR  =  Dun(3) 


INPT  220 
INPT  229 
INPT  230 
INPT  231 


Figure  36.  (continued) 
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XST1( J)=XSTIU(J)  1NPT  521 

rSTI<J)*rS7I0<J)  INPT  522 

3V1  2ST1(J)*ZST1Q< J)  IRPT  523 

G0T02  INPT  524 

t  ***************  ******************************  ********  ************************** 

l***  THIS  PAST  OF  THE  SUBROUTINE  LOADS  THE  FORCES  AND  LOCATIONS  DUE  ****** 

C***  TO  lflPACTS  WITH  TIHE  INTO  THE  'SUMT  MATRIX,  WHICH  IS  CAPAILE  ****** 

C***  OF  50  DISCRETE  TINE  INPUT  VALUES.  THE  MATRIX  SIZE  AS  DEFINED  ****** 

C***  IN  THE  COMMON/DARIER/  ILOCK  IS  (7,50).  THE  SEVEN  POSITIONS  AS  ****** 

L***  THET  ARE  USED  IN  THIS  A0D1FICAT10N  ARE  AS  FOLLOWS:  ****** 

C***  1.  THIS  THE  TINE  LOCATION  ****** 

C*«*  2.  THIS  IS  THE  IMPACT  FORCE  ALIGNED  IN  THE  LOCAL  VEHICLE  ****** 

l**«  X-AXIS  DIRECTION,  POSITIVE  OR  NEGATIVE  ****** 

C**«  3.  THIS  IS  THE  IHPACT  FORCE  ALIGNED  IN  THE  LOCAL  VEHICLE  ****** 

C***  T-AX1S  DIRECTION,  POSITIVE  OR  NEGATIVE  ****** 

C**«  4.  THIS  IS  THE  INPACT  FORCE  ALIGNED  IN  THE  LOCAL  VEHICLE  ****** 

t*«*  Z-AX1S  DIRECTION,  POSITIVE  OR  NEGATIVE  ****** 

C**»  5.  THIS  IS  THE  X  LUCATION  OF  THE  FORCESTRDR  THE  MAES' CENTER-  ****** 

L**«  OF  THE  VEHICLE,  POSITIVE  OR  NEGATIVE  ****** 

C***  6.  THIS  IS  THE  T  LOCATION  OF  THE  FORCES  FROM  THE  NASS  CENTER  ****** 

L*«*  OF  THE  VEHICLE,  POSITIVE  OR  NEGATIVE  ****** 

L***  7.  THIS  IS  THE  2  LOCATION  OF  THE  FORCES  FROM  THE  MASS  CENTER  ****** 

L***  OF  THE  VEHICLE,  POSITIVE  OR  NEBATlVE  ****** 

L**«  JJ  IS  JUST  A  CUUNTER  TO  LOAD  THE  MATRIX  ****** 


l******************************************************************************* 

53  JJ*JJ+1 

suDsn ,  jj)*uuh< i  • 

SUDI ( 2 , J J ) 3OOM ( 2 ) 

SUHK3,  JJ)S0UM(3) 

SUDB<4, JJJ*DUA<4) 

SUDKS,  JJ)*DUM<5) 

SUDD(6, J J)SDUM( A > 

SUBK7,  JJ)=UUM(7 ) 

GO  TO  2 

END  I APT  525 

L  SI  HOLE  VEHICLE  ACCIDENT  SIMULATION  WITH  CURD  IMPACT  -  CRBIMP  CRHP  -1 


Figure  36.  (continued) 


TOO 


NSLCE  =  0 
NUNL0-O 
NUNL 02*0 
YB1VF  =  0.0 


SPUR  Ml 
SFOR  112 
SFOR  113 
SFOK  114 


i;******************************************************************************* 

C«**  THIS  IS  PART  OF  THE  S1QH  IMPACT  MODIFICATION  AS  HADE  BY  JOHN  BETZ  ****** 

C**«  TO  THE  EXISTING  HVQSM  PROORAH.  THE  FOLLOWING  CONAANO  WILL  HAKE 

L**«  THE  PROGRAH  FLOW  SKIP  ALL  THE  PROVISIONS  OF  THIS  PRU6RAA  AND 
C***  DIRECTLY  FIND  THE  IMPACT  FORCES  AND  NOAENTS  AS  FOUND  IN  THE 
C***  MATRIX  'SUBB'.  THE  LINEAR  INTERPOLATION  AETHGD  WILL  BE  USED 
C***  LATER  IN  THIS  SUBROUTINE  TO  FIND  THOSE  FORCES  AND  NOAENTS  AT  ANY 

f***  TINE  INTERVAL.  THE  FORCES  AND  MOMENTS  ARE  LOADED  INTO  THE 

l-***  PROGRAM  NAMES  OF  'SFXS'  ,'SFYS'  ,'SFZS' ,  AND  'SNPS'  ,'SNTS'  ,'SNPSS', 
l***  RESPECTIVELY.  THAT  IS  IN  THE  ORDER  OF  X,  Y,  AND  l. 


****** 

****** 

****** 

****** 

****** 

****** 

****** 

****** 


1FTC0LL.E0.1.)  60  TU  1002 
IF ( INDB.EQ.O)  RETURN 
IB  «  (INOB+11/2 
2  DO  3  1=1,3 


SPUR  115 
SFOR  116 
SPUR  117 


1000  FDKHAT(F7.4, 2F7.T',5T3, T2F8.13' -  - * - SFOR  470 

NLDCTR  =  NLDCTR+1  SFOR  471 

RETURN  SFOR  472 

i;******************  *****************•***********».«****•**•****••*•**•**•*****•** 

C**«  THIS  is  where  the  program  goes  to  USE  THE  impact  forces  anb  ****** 

L***  LOCATIONS  TO  FIND  THE  FORCES  AND  MOMENTS  TO  USE  FOR  THE  CAR  MOTION  ****** 

C***  EQUATIONS.  THIS  USES  LINEAR  INTERPOLATION  AND  OUTPOTS  THE  FORCES  ****** 

C***  AND  MOMENTS  TO  USE  IN  THE  NATRIX  SUBROUTINE.  ****** 

^***t **************  ************************************************************* 

1002  1=2 

1003  IF(T.L£.SUBB< 1 ,1 ) >  GO  TO  1004 
1=1+1 

GO  TO  1903 

1004  IM1*I-1 

VARTMT-SUBBl  1 ,1N1  ;J/<SUBD<  1 ,  It-SUBKl  ,IM1  )> 

SFXS=SUBH(2, !H1 )+VART*(SUBB<2, 1 1-SUBB(2, INI ) > 

SFYS=SUBB(3,IM1 )+VAHT*TSUBB<3,I)*SUBB<3,IM1 > ) 

SFZS=SUBI(4, IMl )+VART*<SUBB<4,ll'SUBB(4,lNl H 
ECKS*SUBB(3,IM1 )+VART*<SUBB(5,l)*SU»B(5rIMl)) 

WHY®  SUBB<6, 1M1 )+VART*(SUBB< 4 , I 1-SUBK6,  INI ) ) 

ZEE®  SUBB <7, IMl )+VAKT*<SUBB<7, 1 )-SUBI(7,lM1 ) ) 

SNPS=SFZS*WHY-SFYS*ZE£ 

SNTS=SFXS*ZEE-SFZS*£CKS 

SNPSS=SFYS*ECKS-SFXS*UHY 

RETURN 

END  SFOR  473 

l*  VGRA  -1 

C  SINGLE  VEHICLE  ACCIDENT  SIMULATION  WITH  CURB  INPACT  -  VGORNT  WORN  0 


Figure  37.  Modification  of  the  SFORCE  Subroutine. 


»«1TE< 4,2028)  APF.APR  1BOT  340 

2028  FURHAT ( 4VH0ANT I-P1TCH  TABLES  FOR  CIRCUMFERENTIAL  TIRE  FURCE  /,  LOOT  341 
1  5H  APF  =  ,  11F10.3/  ,  5X,10F10.3/,5H  APR*,  11F10.3/,5X,10F10.3)  1DUT  342 


CM*  THIS  IS  PART  Of  THE  SIGN  IAPACT  MJBIFICATION  TO  THE  EXISTING  M**** 
C*«*  HVOSH  PRUGRAH  AS  HADE  8Y  JOHN  BETZ.  THE  FOLLOU1N6  COMMANDS  •**»■** 
CM*  PRINT  OOT  THE  IHP ACT  FORCES  A NO  LOCATIONS  AS  READ  IN  IN  THE  ****** 
C***  INPUT  SUBROUTINE.  THE  FORCES  ARE  IN  POUNDS,  THE  DISTANCES  ARE  ****** 
C«*  IN  INCHES,  AND  THE  TIRE  IS  IN  SECONDS.  THIS  ALSO  OSES  THE  'COLL'  ****** 
C**«  TRIGGER  TO  INDICATE  WHETHER  THE  1HPAC7  PART  OF  THE  PROGRAM  IS  ****** 
C***  if»  EFFECT.  ****** 


C** 


***************** *********************************** ************************ 


IF (COLL.NE .  1 . )  GO  TO  1020 
PRINT  1030 

1030  FORMAT!'-', 'SIGN  IMPACT  FORCE  ARD  LOCATION  DATA  *•*******•*•*•• 

1 **************4*****************«****************************' ) 


PRINT  1033 

1033  FORMAT ('O', 'TIME  (S£C)',5X,'X  FORCEILB) ' ,5X, 'Y  FORCE(LB) ' , 

i3x,'z  force(lf)~,7x,'x  <iRr-,?xf'Y-nHr;?Trr  rorrr 

00  1037  I-!,NUMT 

PRINT  1036, SUBS ( 1,1) ,SUBB(2, 1 ) , SUBB(3, I ) , SUIB( 4, I ) ,SUBD(3, I ) 

1 ,SUBB<6,1) ,SUBB(7,I) 

1036  FOKNAT ( '  ',FI0. 3, 6X,F 10. 3,4X ,F10.3,6X,F1 0.3, 5X,F1 0 .3,5X,F1 0.3 , 
15X.F10.3) 

1037  CUNT1NUE 
1020  CONTINUE 

CALL  DRIVID 
WRITE  (6,1007) 

1007  FORMAT  ( 1H1 ) 

RETURN 


IDOT  343 
IDU1  344 
I DO!  345 
IDUf  346 


Figure  38.  Modification  of  the  IDOUT  Subroutine 
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